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ABSTRACT 


Twenty-one  progenies  of  smooth  bromegrass  ( Bromus  inermis  Leyss.)  from  a 
7x7  half  diallel  cross  with  their  parents  were  evaluated  at  four  locations  in  Alberta  and 
in  two  years  for  the  genetic  variation  and  stability  in  expression  of  their  yield  potentials 
along  with  leaf  and  tiller  characteristics. 

Years,  locations  and  their  interactions  were  highly  significant  in  the  combined 
analysis  of  variance,  and  each  of  these  effects  was  therefore  treated  as  an  individual 
environment  in  subsequent  analysis.  The  partition  of  the  genotypic  variance  into 
general  combining  ability  (GCA)  and  specific  combining  ability  (SCA)  showed  that  for 
all  characters  except  spring,  fall  and  annual  yield,  GCA  was  more  important  than 
SCA,  indicating  the  importance  of  additive  genetic  effect.  Although  the  genotype  x 
environment  (GE)  interaction  were  highly  significant,  variation  accounted  for  by 
combining  ability  effects  was  generally  higher  than  the  interaction  effects  of  GCA  and 
SCA  with  environments. 

The  joint  regression  analysis  further  partitioned  the  GE  interaction  into 
heterogeneity  among  regressions  and  its  residual.  A  significant  part  of  the  interaction 
was  ascribed  to  the  heterogeneity  among  regression  lines  for  all  characters  and 
therefore  was  predictable.  However,  the  residual  components  were  significant  in  most 
cases.  Hence  there  was  also  the  presence  of  some  unpredictable  and  unaccountable 
variation. 

Regression  coefficient  and  deviation  from  regression  line  for  each  genotype 
were  the  two  stability  parameters  to  be  considered  together  with  mean  performance  in 
the  evaluation  of  each  genotype.  Five  high  yielding  genotypes,  namely  12,  13,  16,  25 
and  36,  had  general  adaptability,  while  genotypes  23  and  26  were  specially  suited  to 
poor  environment. 

Combining  ability  analysis  revealed  that  GCA  and  SCA  were  both  important 
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in  the  expression  of  mean  yields  (spring,  fall  and  annual).  Inheritance  of  linear 
regression  was  controlled  predominantly  by  GCA  whereas  both  GCA  and  SCA  were 
equally  important  in  the  expression  of  deviation.  Genotype  1  was  the  most  desirable 
parent,  as  it  transmitted  to  its  progeny  high  yield  potential  and  average  linear  response. 
The  other  good  combiner  for  yield  was  the  genotype  6  but  it  transmitted  below  average 
linear  response. 

SCA  analysis  revealed  that  genotypes  12,  13  and  16  combined  both  high  yield 
and  average  stability.  Genotype  34  was  found  superior  in  favourable  environments,  but 
its  non-linear  fluctuations  were  rather  high. 

Phenotypic  and  genotypic  correlation  coefficients  indicated  that  all  the  yield 
types,  including  spring,  fall,  annual  yields,  and  yield  per  area,  were  all  significantly 
correlated.  Very  few  morphological  characters  seemed  to  associate  with  yields,  among 
these  were  tiller  density  and  plant  height.  Beside  these  two  characters,  other  characters 
such  as  leaf  area,  leaf,  stem  and  tiller  dry  weight  were  also  correlated  with  yield  per 
area  but  not  with  other  yield  types. 

Both  direct  and  transformed  scales  were  used  throughout  this  study.  Tests  of 
significance  among  means  were  generally  similar  on  both  scales,  but  for  heterogeneity 
among  regressions,  the  tests  produced  different  results. 

The  results  indicated  that  a  more  complicated  approach,  such  as  recurrent 
selections  involving  multi-location  and  multi-year  testing,  seems  necessary  in  breeding 
for  high  yielding  bromegrass  cultivars  in  Alberta. 
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I.  INTRODUCTION 


Valid  interpretations  of  mechanisms  of  inheritance  as  well  as  predictions  of 
performance  in  breeding  programs  depend  mainly  on  accurate  assessments  of  genotype 
values.  These  assessments  must  be  made  from  data  on  phenotypes  that  reflect  both 
nongenetic  and  genetic  influences  on  plant  development  (Moll  and  Stuber,  1974; 
Comstock  and  Moll,  1963).  Unfortunately  the  genetic  effects  are  not  independent  of  the 
nongenetic  environmental  effects.  The  phenotypic  response  to  a  change  in  environment 
is  not  the  same  for  all  genotypes;  these  are  reflected  by  the  difference  in  the  relative 
rankings  of  genotypes  in  different  environments.  This  interplay  of  genetic  and 
nongenetic  effects,  i.e.  of  genotype-environment  (GE)  interaction,  reduces  the 
correlation  between  genotype  and  phenotype,  which  in  turn  reduces  confidence  in 
inferences  from  experimental  data  relevant  to  both  plant  improvement  and  inheritance 
mechanisms. 

According  to  Allard  and  Bradshaw  (1964),  environment  variation  can  be 
classified  into  two  types:  predictable  and  unpredictable.  Predictable  variation  includes 
all  permanent  characters  of  environment,  such  as  climate  and  soil  type,  as  well  as  cyclic 
tluctuations  such  as  day  length.  Also  included  in  this  category  are  those  which  can  be 
fixed  at  will  (such  as  planting  date,  sowing  density,  harvest  methods,  and  other 
agronomic  practices).  Unpredictable  variation  includes  fluctuations  in  weather,  such  as 
amount  and  distribution  of  rainfall  and  temperatures.  Although  the  distinction  between 
the  two  categories  is  not  always  clear  cut,  they  have  distinctly  different  impacts  in 
breeding  programs,  both  on  operational  procedures  at  the  selection  stages  and  on 
testing  stages.  Therefore,  results  based  on  one  location  and  one  year  have  limited 
practical  implication. 

The  significance  of  GE  interactions  to  plant  breeders  depends  on  their 
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objectives.  If  cultivars  that  perform  well  over  a  broad  spectrum  of  environments  are 
desired,  small  GE  interactions  and  well-buffered  cultivars  are  favoured.  If  cultivars  that 
are  adapted  to  very  specific  or  predictable  environments  are  desired,  large  interactions 
may  be  beneficial,  whereas  buffering  may  be  of  little  consequence.  There  are,  however, 
dangers  in  an  extremely  inflexible  cultivar  because  of  both  the  state  of  influx  with 
regards  to  management  ideals  and  the  mismanagement  which  can  be  expected  to  occur 
on  some  occasions  under  practical  farming  conditions  and  the  variable  climate. 
Information  on  performance  in  a  range  of  environments  is  of  vital  importance  to  the 
breeder,  whether  his  interest  is  in  the  production  of  a  cultivar  with  a  narrow  or  with  a 
relatively  wide  spectrum  of  adaptation.  In  Alberta,  the  varying  soil  climatic  patterns 
from  the  north  to  the  south,  and  from  the  west  to  the  east  are  so  great  that  perhaps 
large  interactions  and  cultivars  adapted  to  more  specific  regions  of  the  province  may  be 
more  feasible  than  well-buffered  cultivars  for  the  whole  province. 

The  ability  of  individuals  to  respond  to  varying  environments  by  changing 
phenotypic  expression  is  of  special  interest  in  herbage  plants,  since  no  other  crops  are 
expected  to  yield  under  such  a  diversity  of  climatic,  edaphic  and  management 
conditions.  This  is  especially  true  in  perennial  grasses,  where  successive  harvests  are 
subject  to  marked  contrasts  in  the  external  environment.  There  are  also  certain  areas  in 
Canada  where  high  protein  legumes  cannot  be  grown  (Wilson  and  Winch,  1977).  Under 
such  conditions  grass  will  remain  the  prime  component  of  programs  for  pasture  or 
stored  feed,  but  adequate  amounts  of  nitrogen  must  be  applied  to  obtain  high  yields. 
Cultivars  or  genotypes  which  have  a  certain  amount  of  stability  but  will  be  responsive 
to  favourable  environment  (e.g.  nitrogen)  would  be  favoured.  Thus,  GE  interaction 
studies  are  of  the  utmost  practical  importance  in  measuring  the  genotypic  reaction  to 
environmental  variations. 

Until  recently,  differential  responses  of  genotypes  or  varieties  to  varying 
environments  were  determined  using  the  analysis  of  variance  (Comstock  and  Robinson, 
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1952).  Over  the  past  few  years  interest  in  the  form  of  regression  analysis  first  proposed 
by  Yates  and  Cochran  (1938)  has  been  revived  by  the  studies  of  Finlay  and  Wilkinson 
(1963)  and  the  Birmingham  School  of  Genetics  (e.g.  Bucio  Alanis  and  Hill,  1966; 
Perkins  and  Jinks,  1968).  The  significant  feature  in  the  regression  analysis  is  that 
genetic  expression  can  be  simply  and  predictably  related  to  the  environment  when  the 
latter  is  measured  by  its  effect  on  the  character  under  study. 

The  present  investigation  was  initiated  to  examine  the  magnitude  of  GE 
interaction  with  a  7x7  half  diallel  cross  of  smooth  bromegrass  ( Bromus  inermis  Leyss.), 
which  had  been  studied  for  a  number  of  aspects,  including  leaf  and  tiller  characteristics 
(Tan  et  al.,  1976a,  b;  Walton,  1976;  Murchison,  1977),  canopy  (Tan  et  ah,  1977),  and 
quality  (Tan  et  ah,  1978).  More  precisely,  the  objectives  are: 

1.  to  evaluate  the  magnitudes  of  general  combining  ability,  specific  combining  ability 
and  their  interactions  with  different  environments; 

2.  to  use  the  regression  procedure  to  analyse  and  study  the  GE  interaction  on  the 
data  obtained.  The  estimates  of  linear  and  non-linear  parameters  provide  an 
account  of  the  dynamic  response  of  bromegrass  genotypes  to  changing 
environments,  and  may  be  used  with  mean  performance  to  assess  the  potentialities 
of  various  genotypes;  and 

3.  to  relate  stability  parameters  with  combining  ability  effects  in  the  selection  of 
bromegrass  genotypes  which  combine  high  yield  with  average  stability  in  differing 


Alberta  environments. 


II.  LITERATURE  REVIEW 


Early  Approaches  to  the  Study  of  GE  Interactions  -  the  Use  of  Variance  Components 

The  existence  of  interactions  between  genotypes  and  environmental  factors  has 
long  been  recognised.  Factorial  design  was  adapted  to  the  analysis  of  GE  interaction, 
because  the  total  variation  ascribable  to  genotypes  and  environments  was  partitioned 
into  three  orthogonal  comparisons:  (1)  measuring  the  differences  between  genotypes,  (2) 
measuring  environmental  differences,  and  (3)  assessing  their  joint  effects  (Hill,  1975). 

Fisher  laid  the  foundations  for  use  of  factorial  design  and  analysis  of  field 
experiments  (Hill,  1975).  Among  those  who  developed  it  further  and  were  responsible 
for  using  the  analysis  of  variance  to  investigate  GE  interactions  are,  e.g.  Sprague  and 
Federer  (1951),  Comstock  and  Robinson  (1952),  Hanson,  Robinson  and  Comstock 
(1956)  and  Comstock  and  Moll  (1963).  Unbiased  estimates  of  the  genetic  and  GE 
components  of  variance  can  be  readily  obtained  by  equating  the  expected  mean  squares 
with  those  calculated  from  the  experiment.  The  precise  form  taken  by  the  expected 
values  of  the  mean  squares  will  depend  upon  the  underlying  assumptions  made  in  the 
analysis.  Unbiased  estimates  of  the  components  of  variance  may  be  attained  through 
properly  designed  experiments.  Accuracy  of  these  estimates  will  depend  on  the  size  of 
the  experiments. 

The  Use  of  Regression  Methods 

The  idea  of  breaking  down  an  interaction  into  several  parts  was  missing  in  the 
variance  component  approach.  The  method  of  partitioning  GE  interaction  was  first 
devised  by  Yates  and  Cochran  (1938),  but  attracted  little  attention  until  Finlay  and 
Wilkinson  (1963)  rediscovered  the  same  method  and  used  it  for  an  analysis  of 
adaptation  in  a  trial  with  277  varieties  of  barley  in  seven  environments.  A  similar 
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method  was  also  used  by  Walton  in  the  study  of  cotton  (1957;  1958;  and  1961).  This 
method  involves  computing  for  each  variety  the  regression  of  individual  yield  on  the 
environmental  mean,  which  was  measured  by  the  mean  yield  of  all  varieties  for  each 
site  and  season.  The  regression  has  since  emerged  as  one  of  the  most  useful  statistical 
tools  for  analysing  GE  interactions. 

The  approach  falls  into  two  parts,  a  conventional  analysis  of  variance  being 
followed  by  a  joint  regression  analysis.  From  the  joint  regression  analysis,  the  GE 
interaction  sum  of  squares  can  be  partitioned  into  two  orthogonal  items,  heterogeneity 
of  regression  and  deviation  from  regression.  The  former  measures  that  portion  of  the 
GE  interactions  which  is  due  to  differences  between  the  fitted  regression  lines,  with  (t-1) 
degrees  of  freedom  (d.f.)  (t  being  the  number  of  genotypes);  the  latter  measures  the 
accumulated  deviations  of  the  observed  values  around  these  fitted  lines,  with 
(t-l)(s-2)d.f.  (s  being  the  number  of  environments).  Each  of  these  terms  can  be 
compared  with  the  experimental  error,  and  the  heterogeneity  of  regression  further 
compared  with  the  deviation  in  order  to  see  if  it  accounts  for  a  significantly  large  part 
of  the  observed  interaction. 

Regression  methods  were  also  considered  by  Rowe  and  Andrew  (1964)  and 
Eberhart  and  Russell  (1966),  who  added  together  the  sums  of  squares  for  environments 
and  GE  interactions  and  repartitioned  this  into  a  linear  component  between 
environments  with  1  d.f.,  a  linear  component  of  the  GE  interaction  with  (t-l)d.f.,  and 
deviations  from  regression.  The  deviations  were  found  separately  for  each  of  the  t 
genotypes  with  (s-2)d.f.  each. 

The  linear  regression  technique  is  versatile  and  can  be  used  to  predict  the 
performance  either  of  genotypes  in  environments  other  than  those  sampled 
experimentally  or  of  segregating  generations  from  the  non-segregating  generations  from 
which  they  were  derived.  The  technique  has  been  usefully  applied  to  a  number  of 
different  crop  and  plant  species  since  its  rediscovery,  such  as  barley  (  Hordeum  vulgare 
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L.)  (Finlay  and  Wilkinson,  1963;  Paroda  and  Hayes,  1971),  wheat  (  Triticum  aestivum 
L.)  (Kaltsikes  and  Larter,  1970),  maize  (  Zea  mays  L.)  (Eberhart  and  Russell,  1966; 
Ottaviano  and  Sari  Gorla,  1972;  Dhillon  and  Singh,  1977),  orchard  grass  (  Dactylis 
glomerata  L.)  (Breese,  1969);  perennial  ryegrass  (  Lolium  perenne  L.)  (Troughton,  1970; 
Hill  and  Samuel,  1971),  soybean  (  Glycine  max  L.  Merr.)(Baihaki,  Stucker  and 
Lambert,  1976),  strawberry  (  Fragaria  chiloensis )  (Gooding,  Jennings  and  Topham, 
1975),  peas  (  Pisum  sativum  L.)  (Snoad  and  Arthur,  1974),  Nicotiana  rustica  (Bucio 
Alanis,  1966;  Bucio  ALanis  and  Hill,  1966;  Perkins  and  Jinks,  1968),  Arabidopsis 
thaliana  (Westerman  and  Lawrence,  1970),  and  Schizophy/lum  commune  (Fripp  and 
Caten,  1971).  In  many  instances,  the  linear  regression  technique  adequately  describes 
the  behaviour  of  genotypes  over  a  range  of  environments. 

The  linear  regression  technique  has  been  widely  adapted  in  various  designs  and 
experiments  in  recent  years.  For  example,  Breese  and  Hill  (1973),  Hill  (1973),  and 
Wright  (1971,  1977)  have  adapted  the  techinque  to  the  analysis  of  competition 
experiments  based  upon  the  diallel  arrangement,  that  is  where  the  competitors  are 
grown  as  monocultures  in  all  possible  50:50  mixtures. 

Lin,  Binns  and  Thompson  (1977)  described  a  model  combining  the  features  of 
Griffing's  diallel  cross  analysis  with  regression  analysis  for  genotype-environment 
interactions.  They  claimed  the  main  advantage  for  studying  GE  interactions  in  a  diallel 
cross  experiment  is  that  ‘both  the  combining  abilities  of  genetic  effects  and  their  linear 
function  of  combining  ability  components  by  environment  interaction  can  be  studied 
simultaneously’.  Furthermore,  the  authors  (Lin  et  al.,  1977)  claimed  the  model  provides 
not  only  a  direct  and  easy  biological  interpretation  but  also  an  easy  assessment  of  both 
heterosis  and  heterosis  by  environment  interaction. 
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Limitations  of  Regression  Methods  on  the  Use  of  Environmental  Mean 

Although  the  linear  regression  approach  has  been  very  widely  used,  it  has 
frequently  been  the  subject  of  controversy.  Freeman  and  Perkins  (1971)  have  criticized 
the  use  of  marginal  means  of  the  environments  as  independent  variables  in  the 
regression  analysis,  claiming  their  use  violates  the  fundamental  statistical  assumptions 
of  regression  analysis.  If,  as  Freeman  and  Perkins  (1971)  further  suggested,  it  is 
necessary  that  points  on  the  abscissa  (  x-variate  )  should  be  known  precisely  and  occur 
at  equal  intervals  over  the  observed  range,  such  as  controlled  environment  chambers, 
then  the  results  would  inevitably  be  of  more  theoretical  than  practical  interest. 

However,  Hardwick  and  Wood  (1972)  proposed  that  if  a  large  number  of  genotypes  are 
included,  and  if  the  environmental  range  is  such  that  the  between-environments  mean 
squares  are  significantly  greater  than  the  error  mean  square,  any  bias  which  results 
should  not  prove  serious  in  practice. 

In  view  of  the  criticism  of  Freeman  and  Perkins  (1971),  different  ways  of 
providing  an  independent  assessment  of  the  environment  have  been  suggested: 

1.  The  use  of  control  genotypes  or  extra  replicates  of  the  full  genotype  set  (Fripp  and 
Caten,  1971;  Jinks  and  Connolly,  1973).  Fripp  and  Caten  (1971)  found  that  with  a 
large  number  of  genotypes,  the  results  for  regression  on  a  control  genotype  differed 
little  from  those  for  regression  on  the  mean  of  the  test  genotypes.  Further,  Fripp 
(1972)  compared  both  biological  and  physical  measures  of  the  environment  and 
found  that  the  analyses,  again  for  a  large  number  of  genotypes,  gave  very  similar 
results  for  all  reasonable  external  measures  and  the  environmental  mean. 

2.  The  use  of  parental  genotypes  as  standards  in  relation  to  any  generation  derived 
from  crosses  between  them  (Bucio  Alanis  and  Hill,  1966;  Breese,  1969;  Bucio 
Alanis,  Perkins  and  Jinks,  1969;  Jinks  and  Perkins,  1970;  Perkins,  1970).  Perkins 
and  Jinks  (1973)  regressed  values  for  a  large  number  of  genotypes  on  their 
environmental  means  and  on  values  derived  from  other,  closely  related,  sets  of 
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genotypes;  they  found  that  all  analyses  generally  gave  similar  values  for 
significance,  but  that  regressions  on  means  derived  from  only  a  few  independent 
genotypes  were  sometimes  so  insensitive  as  to  give  rise  to  problems  of 
interpretation. 

3.  By  regressing  the  performance  of  the  i th  genotype  onto  an  index  composed  of  the 
remaining  genotypes  (Mather  and  Caligari,  1974).  This  removes  the  statistical 
objection  without  needing  to  include  extra  individuals  in  the  experiment,  though 
estimates  so  obtained  will  be  distorted,  both  by  error  variation  and  by  any 
departure  from  linearity  on  the  part  of  the  individual  regressions. 

Finally,  Tai(  1971)  has  employed  'structural  relationship  analysis'  to  overcome 
the  limitations  of  regressing  one  set  of  variables  on  another  which  is  not  independent 
of  it.  Assuming  that  the  environmental  and  GE  effects  are  jointly  normally  distributed, 
it  is  possible  to  derive  maximum  likelihood  estimates  of  the  linear  response  of  the  /th 
line  to  the  environmental  effects  (a,)  and  the  deviation  from  the  linear  response  in 
terms  of  the  magnitude  of  the  error  variance)  AT  These  estimates  may  be  compared 
with  their  counterparts  obtained  by  regression  analysis. 

Stability  Parameters 

Considerable  interest  exists  in  the  mechanisms  by  which  an  individual  stabilizes 
its  behaviour  in  the  face  of  varying  environmental  influences  (Bradshaw,  1965).  The 
definitions  of  stability  are  many  and  varied,  even  within  the  confines  of  GE 
interactions. 

Plaisted  and  Peterson  (1959)  based  their  measures  of  stability  upon  the 
contribution  of  the  /th  genotype  to  the  GE  interaction  sum  of  squares.  Finlay  and 
Wilkinson  (1963)  defined  stable  genotypes  as  those  with  a  low  regression  coefficient 
(6,<1.0)  while  genotypes  with  a  high  regression  value  (£>,>1.0)  are  unstable.  Stable 
varieties  with  high  mean  yields  are  regarded  as  those  specifically  adapted  to 
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low-yielding  environments,  and  unstable  varieties  with  high  mean  yields  are  specifically 
adapted  to  high-yielding  environments.  Hanson  (1970)  devised  a  composite  measure  of 
stability  which  combines  the  contribution  of  the  / th  genotype  to  the  GE  interaction  sum 
of  squares  with  its  response  to  environmental  change.  Eberhart  and  Russell  (1966)  also 
used  bi  as  the  first  measure  of  stability  but  went  further  and  regarded  the  mean  square 
of  deviation  Sd2  as  a  second  measure.  They  defined  a  stable  genotype  as  a  line  which 
has  Z?j=  1 .0  and  Sd2—  0.  Perkins  and  Jinks  (1968)  obtained  similar  parameters  to  that  of 
Eberhart  and  Russell  (1966)  and  Finlay  and  Wilkinson  (1963)  from  a  biometrical 
genetic  model,  but  their  regression  coefficients  /3,( i.e.  6,-1)  centered  around  0.0  rather 
than  1.0. 

Tai  (1971)  used  a,  as  one  measure  of  stability  and  also  defined  a  second 
measure  A,.  Tai's  parameters,  cti  and  A,,  have  equivalent  meaning  to  the  6i  and  Sd,2,  of 
Eberhart  and  Russell  (1966),  respectively,  except  that  Tai's  a;  and  Ai  center  around  0.0 
and  1.0,  respectively.  Tai  (1971)  defined  ‘a  perfectly  stable  variety  was  (a,, A,)  =  (-1,1) 
and  a  variety  with  average  stability  was  («,>,)  =  (0,1)’. 

Scales  and  Transformations 

Most  estimates  of  stability  parameters  have  been  made  on  a  direct  scale,  but 
square  root  (Westerman  and  Lawrence  1970)  and  logarithmic  scales  (Finlay  and 
Wilkinson,  1963)  have  been  also  used.  Jowett  (1972)  estimated  stability  parameters  for 
grain  sorghum  ( Sorghum  bicolor  L.)  on  both  direct  and  logarithmic  scales.  The  ranking 
of  single  crosses,  three-way  crosses,  and  varieties  by  regression  coefficients  differed  on 
the  two  scales.  This  supports  the  observation  by  Knight  (1970),  that  different 
conclusions  may  be  drawn  from  regression  coefficients  when  measured  on  different 
scales.  Breese  and  Hill  (1973)  also  found  that,  not  only  did  a  log  transformation  of 
data  taken  from  competition  experiments  fail  to  remove  the  interactions  which  were 
present,  but  it  also  completely  altered  the  interpretation  of  the  results.  By  contrast. 
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Mather  (1971)  asserted  that  no  deep-seated  biological  significance  can  be  attached  to 
the  particular  scale  chosen. 

Recent  evidence  on  oats  ( Avena  sativa  L.)  by  Eagles,  Hinz  and  Frey  (1977) 
further  supports  the  postulation  that  regression  coefficient  parameters  differed 
substantially  on  the  two  scales  i.e.  direct  and  transformed  scales.  Therefore,  they 
warned  that  great  care  should  be  taken  when  comparing  regression  coefficient 
parameters  estimated  on  different  scales.  However,  they  added  that  if  the  GE 
interaction  can  be  removed  by  using  any  of  the  family  of  transformations  Y=Xa  where 
a  ranged  from  0  to  1  (or  any  other  monotonic  transformation),  the  interaction  will 
likely  be  of  the  type  which  cannot  be  exploited  by  selection. 

The  Use  of  Multivariate  Techniques  to  Study  GE  Interactions 

There  can  be  no  doubt  that  the  most  useful  of  the  GE  interaction  techniques 
has  been  the  regression  approach.  For  successful  application  of  the  regression  method, 
a  very  high  proportion  of  the  interaction  sum  of  squares  should  be  explained  by  linear 
regression.  The  conditions  making  for  success,  i.e.  linearity  of  regression,  are  very 
difficult  to  determine,  and  one  set  of  characters  has  frequently  been  found  to  give 
linear  regressions,  while  another  set  of  characters  measured  on  the  same  set  of 
genotypes  has  not  (Freeman,  1973). 

In  practice  there  are  often  wide  deviations  from  linearity.  There  are  situations 
where  the  linear  regression  technique  may  over-simplify  the  true  response  pattern  to  an 
extent  which  could  lead  to  erroneous  conclusions  (Witcombe  and  Whittington,  1971). 
The  results  of  a  number  of  studies  suggest  that  there  is  often  more  than  one  way  in 
which  responses  differ,  that  is,  the  interaction  contains  more  than  one  significant 
principal  component.  Multivariate  methods  will  therefore  likely  be  alternative 
techniques  to  the  analyses  of  GE  interactions. 

Multivariate  techniques  of  anlayses  are  essentially  an  extension  of  the 
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univariate  techniques.  The  purposes  of  multivariate  analysis  may  be  summarized  in 
terms  of  the  analysis  of  GE  interactions  as  follows  (Hill,  1973): 
first,  to  assess  the  simultaneous  effects  of  a  number  of  environmental  factors  where 
these  can  be  measured  and  ranked  in  order  of  importance  by  determining  how  much  of 
the  observed  variation  is  accounted  for  by  each  individual  factor,  or  composite  factor 
derived  therefrom; 

secondly,  to  maximize  differences  between  cultivars  (or  environments)  relative  to 
differences  within  cultivars  (or  environments). 

Multivariate  techniques  have  not  been  widely  used  in  plant  breeding,  still  less 
in  the  analysis  of  GE  interactions.  Freeman(  1973)  listed  several  multivariate  techniques; 
amongst  them  principal  component  analysis,  canonical  analysis,  cluster  analysis,  and 
factor  analysis,  which  could  be  of  prospective  value  in  the  analysis  of  GE  interactions. 
Tai  (1975)  also  proposed  the  method  of  path  coefficient  based  on  postulated  causal 
relationships. 

Freeman  and  Dowker  (1973)  applied  principal  component  analysis  to  data 
recorded  from  a  series  of  yield  trials  in  carrots  after  the  joint  regression  analysis  had 
been  only  partially  successful  in  explaining  the  observed  GE  interactions.  They 
concluded  the  use  of  principal  component  analysis  supplied  no  additional  information 
beyond  that  obtained  from  the  analysis  of  variance.  A  similar  conculsion  was  drawn  by 
Perkins  (1972)  using  the  same  analysis.  Shukla  (1972;  cited  by  Freeman,  1973)  also 
demonstrated  no  clear  advantage  of  canonical  analysis  over  regression  analysis  on 
environmental  variables. 

Multivariate  techniques  will  undoubtedly  figure  prominently  in  the  GE 
interaction  context,  particularly  where  linearity  fails.  Even  so  the  linear  regression 
technique  will  continue  to  play  an  important  part  in  furthering  our  understanding  of 
GE  interactions  because,  despite  its  imperfections,  it  does  have  the  twin  merits,  of 
simplicity  and  biological  relevance  (Hill,  1975).  There  are  without  doubt  biological 
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problems  which  can  be  most  effectively  solved  by  multivariate  techniques.  But,  as  Hill 
(1975)  stated,  'there  is  the  very  real  danger  that  biological  relevance  will  be  sacrificed 
for  statistical  pedantry’,  against  which  we  must  guard. 


III.  MATERIALS  AND  METHODS 


Selections  of  Environments  and  Parents,  and  Experimental  Designs 

A  seven-parent  diallel  cross,  without  reciprocals,  of  smooth  bromegrass 
( Bromus  inermis  Leyss.)  was  grown  at  four  locations  over  three  years.  The  locations 
were  chosen  to  provide  differences  in  soil  type,  annual  rainfall,  and  seasonal 
temperatures  representing  the  cultivated  areas  within  the  province  of  Alberta.  The 
regions  and  locations  were:  Northwest  -  Beaverlodge,  Research  Station,  Canada 
Department  of  Agriculture;  Central  -  Edmonton,  Parkland  farm;  Central-east  - 
Kinsella,  University  Ranch;  and  South  -  Lethbridge,  Research  Station,  Canada 
Department  of  Agriculture.  The  soil  characteristics,  fertility  and  climatic  factors  for 
each  site  are  described  in  Tables  1,  2,  and  3.  All  locations  had  been  fallowed  the  year 
before  planting.  Approximately  88  and  68  kg/ha  of  nitrogen,  and  55  and  36  kg/ha  of 
phosphorous  were  applied  in  September  1975  and  1976,  respectively,  to  all  locations. 

The  details  of  the  origin  and  selections  of  the  parental  clones  are  described  in 
Table  4  .  Seeds  of  the  21  single  crosses  were  obtained  by  mutual  pollination  (Mishra 
and  Drolsom,  1972)  without  emasculation  in  the  greenhouse  during  the  winter  of 
1972-73.  Hybrid  seed  and  tillers  of  the  parents  were  planted  in  the  green  house  in 
individual  peat-pots  for  a  month.  The  parental  propagules  and  Fi  progenies  were  then 
transplanted  to  each  site  in  1975  (actual  planting  dates  for  each  site  are  given  in  Table 
5)  in  randomized,  complete  block  design  with  six  replications.  Because  of  the  limited 
number  of  hybrid  seed,  each  plot  consisted  of  five  plants  spaced  120  cm  between  plants 
within  and  between  plots.  The  test  was  bordered  by  rows  of  cultivar  ‘Magna’  with  the 
same  spacing  as  that  of  the  plots.  Cultivation  was  limited  to  between  plots,  so  that 
plants  within  plots  naturally  grew  into  each  other  after  two  years  and  were  considered 
as  a  sward  in  the  third  year  (1977).  The  trials  were  not  harvested  until  1976,  the  second 
year  of  establishment. 
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Table  1.  Description  of  soil  of  each  of  the  four  locations  in  which  the  smooth  bromegrass  genotypes  were 
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on  irrigated  site.  Irrigated  twice  during  the  experimental  period: 
i.  Four  weeks  after  establishment  in  June  1975, 
ii.  Two  weeks  before  1st  harvest  in  June  1977,  equivalent  to  approximately  3  inches  of  water 
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Table  2.  Descriptions  of  initial  and  final  soil  fertility  levels  of 
the  bromegrass  plots  at  four  locations  and  two  dates. 


Depth  of 
Sample 

Location  (cm) 

Fertility  level 

Nitrogen 

(N) 

Phosphorous 

(P) 

Potassium 

(K) 

Soil  Reaction 
(pH) 

I.  September  1975 

kg/ha 

(Initial) : 

Beaverlodge 

0-15 

81 

19 

543 

6.4 

15-30 

53 

3 

400 

5.  8 

Kinsella 

0-15 

58 

34 

511 

6.1 

15-30 

44 

8 

316 

6.5 

Lethbridge 

0-15 

84 

39 

623 

7.9 

15-30 

71 

8 

420 

7.7 

Edmonton 

0-15 

85 

34 

669 

5.9 

15-30 

94 

25 

601 

6.0 

II.  September  1977 

(Final)  : 

Beaverlodge 

0-15 

9 

23 

474 

6.0 

15-30 

4 

3 

477 

5.4 

Kinsella 

0-15 

11 

43 

394 

6.1 

15-30 

2 

8 

260 

6.5 

Lethbridge 

0-15 

60 

48 

680 

7.9 

15-30 

94 

11 

436 

7.8 

Edmonton 

0-15 

9 

51 

1081 

6.1 

15-30 

5 

37 

676 

6.3 

Table  3.  Description  of  the  major  climatic  factors  for  the  four  locations  in  the  year 
of  establishment  (1975)  and  the  two  following  years. 
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Table  4. 


Parent 

UA5 

UA9 

UA10 

UA12 

B40 

B42 


A  description  of  the  origin  of  the  parental  clones  used  in 
the  diallel  cross. 


Code  Origin 

1  Random  selection  from  'Magna",  an  intermediate  northern- 
southern  type  bromegrass. 

2  Random  selection  from  'Carlton',  a  northern- type 
bromegrass. 

3  Random  selection  from  S7388,  a  synthetic  susceptible 
to  Selenophoma  bromigena  and  Pyrenophora  bromi  from 
Saskatoon  by  Smith  and  Knowles  (Walton,  1974). 

4  Random  selection  from  'Lincoln',  a  southern- type 
bromegrass . 

5  Selection  by  Smith  and  Knowles  (Walton,  1974)  showing 
resistance  to  Selenophoma  bromigena  and  Pyrenophora 
bromi.  Yielded  low  in  a  previous  diallel  experiment. 

6  Selection  by  Smith  and  Knowles  (Walton,  1974)  showing 
resistance  to  Selenophoma  b romigena  and  Pyrenophora 

b romi .  Highest  yielding  resistant  strain  from  a 
previous  diallel  experiment  where  it  also  gave  a  good 
GCA  and  high  yielding  progeny  when  crossed  with  UA9. 

7  Collected  from  an  old  auction  yard  on  a  farm  near 
Sedgewick  in  1970.  Outyielaed  'Carlton'  by  19%  in 
1972  and  17%  in  1971. 
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Table  5.  Records  of  dates  of  planting  and  harvesting  for  the  four  locations. 
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Sampling  Techniques 

Ten  headed  tillers,  two  from  each  plant,  were  randomly  sampled  from  each 
plot  in  June  each  year  prior  to  the  spring  harvests  for  each  location.  Number  of  green 
leaves  per  tiller  was  counted  and  the  leaf  blade  area  per  tiller  was  measured  using  an 
automatic  leaf  area  meter  (Hayashi  Denko,  Model  AAM-5).  The  tillers  were  separated 
into  leaves  and  stems  to  obtain  per  tiller  leaf  dry  weight,  stem  dry  weight,  and  whole 
tiller  dry  weight.  Specific  leaf  weight  (SLW)  was  obtained  by  dividing  leaf  dry  weight 
by  leaf  area. 

Both  fertile  and  sterile  tillers  in  an  area  of  196  cm"  (i.e.  14cm  x  14cm)  were 
taken  from  each  plant  for  all  locations  except  Beaverlodge,  where  the  area  was  reduced 
to  100  cm2  (10cm  x  10cm).  The  sample  was  weighed  to  determine  yield  per  area 
(g/dm:),  and  was  counted  for  the  number  of  tillers.  The  records  were  later  converted  to 
a  per  unit  area  basis.  Plant  height  was  recorded  for  each  plant  before  the  spring 
harvest.  Two  hay  crops  were  harvested  before  anthesis  for  each  of  the  two  years,  the 
dates  are  given  in  Table  5.  Five  whole  plants  from  each  plot  were  harvested  in  1976 
whereas  a  plot  size  of  60cm  x  450cm  was  harvested  in  1977.  Forage  yield  was 
expressed  as  dry  weight  in  grams  per  plot. 

The  harvesting  dates  as  listed  in  Table  5  were  kept  as  close  to  one  another  as 
possible.  However,  some  delay  due  to  unfavourable  weather  was  unavoidable. 


20 


Statistical  Analyses 


Transformation 

One  of  the  assumptions  underlying  the  analysis  of  variance  is  that  the 
experimental  factors  shall  have  the  same  variance.  Error  mean  squares  from  the 
variance  analyses  of  experiments  with  identical  treatments,  but  conducted  in 
different  environments,  are  often  heterogeneous  (Yates  and  Cochran,  1938).  When 
experimental  means  and  error  variances  are  correlated,  a  transformation  to  stabilize 
the  error  variance  is  recommended  before  conducting  a  combined  analysis  (Bartlett, 
1947). 

To  remove  the  correlation  between  the  error  mean  squares  and 
experimental  means,  a  transformation  of  the  form  as  outlined  was  made. 

Y  =  g(X) 

where  X  is  the  character  in  the  original  scale,  Y  is  the  character  on  the  transformed 
scale, 

g(X)  =  Xa 

where  a  is  a  constant  (Hinz  and  Eagles,  1976).  The  method  described  by  Hinz  and 
Eagles  (1976)  was  followed  to  obtain  the  value  of  a.  How'ever,  in  cases  where  the 
Hinz  and  Eagles  method  was  unsuccessful,  a  logarithmic  transformation  was  used. 
Table  6  listed  the  types  of  transformation  applied;  data  for  leaf  number  and  height 
were  homogeneous  and  were,  therefore,  not  transformed. 
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Table  6.  A  list  of  the  type  of  transformation  applied  to  the  characters 
studied  and  the  correlation  coefficients  (r)  between  means  and 
error  mean  squares  of  eight  environments  for  each  character. 


Type  of  Transformation 

Root  Logarithm 


_ r _ 

Before  After 

Transformation  Transformation 


Yield 


Spring  yield 

0.  685 

— 

0.49 

0.04 

Fall  yield 

0.5 

— 

0.89 

0.15 

Annual  yield 

0.  5 

— 

0.59 

0.10 

Yield/area 

— 

(-1)  x 

loge 

0.92 

-0.21 

Leaf  Characters 

Leaf  area 

0.445 

— 

0.97 

0.13 

Leaf  no 

— 

— 

-0.  08 

— 

Leaf  dry  wt 

0.25 

— 

0.67 

-0.17 

Specific  leaf  wt 

— 

(-1)  X 

loge 

-0.07 

0.20 

Tiller  Characters 

Stem  dry  wt 

0.25 

— 

0.91 

-0.05 

Tiller  dry  wt 

0.25 

— 

0.  89 

-0. 10 

Tiller  density 

0.58 

— 

0.45 

-0.07 

Leaf  stem  ratio 

— 

(-1)  X 

iog 

-0.19 

0.33 

Height 

— 

— 

6 

0.14 

— 
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Analysis  of  Variance  and  Dial  lei  Analysis 

All  statistics  were  analysed  using  the  plot  means  of  the  characters.  In  the 
combined  analysis  of  variance,  all  sources  including  years,  locations,  and  genotypes 
were  assumed  fixed,  therefore  the  main  effects  and  interactions  were  tested  against 
their  corresponding  error  variance. 

The  source  of  variation  caused  by  genotypes  (including  parents  and  Fi 
progenies)  in  the  analysis  of  variance  was  further  partitioned  into  general 
combining  ability  (GCA)  and  specific  combining  abiltiy  (SCA).  As  highly 
significant  differences  were  obtained  for  years-locations  combination  for  almost 
every  character  (Tables  13,14,15),  each  year-location  combination  was  considered  an 
environment  in  subsequent  analysis.  The  genotypes  x  years  and  genotypes  x 
locations  were  pooled  to  become  genotypes  x  environments,  which  was  partitioned 
into  GCA  x  environments  and  SCA  x  environments  effects.  The  data  were  analysed 
using  the  computer  program  DIALL  (Schaffer  and  Usanis,  1969). 

Griffing's  (1956)  Model  I,  Method  2  (parents  and  one  set  of  F i s)  and 
Method  4  (one  set  of  Fis)  were  used  to  evaluate  GCA  and  SCA  effects  for  each 
individual  environment. 

Phenotypic  and  genotypic  correlation  for  all  possible  combinations  of  the 
characters  measured  were  calculated  according  to  the  method  outlined  by  Johnson, 
Robinson  and  Comstock  (1955).  Correlation  coefficients  for  GCA  and  SCA  were 
also  computed  from  the  variance  and  covariance  components  as  described  by 
Griffing  ( 1956). 

Regression  Analysis 

Regression  technique  was  used  for  further  analysis  of  the  GE  interactions 
using  the  methods  of  Perkins  and  Jinks  (1968),  and  Freeman  and  Perkins  (1971). 
Two  measures  of  environment  were  used.  The  first  measure  was  the  average  of  all 
entries  in  the  trial  at  that  environment.  Each  of  the  28  genotypes  was  regressed 
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upon  that  environmental  mean.  The  second  measure  was  the  average  of  the  seven 
parental  genotypes,  the  analysis  was  performed  by  regressing  the  mean  of  each  of 
the  21  Fi  progenies  on  the  environmental  mean. 

All  factors  were  again  assumed  fixed.  All  sources  including  heterogeneity  of 
regression  and  its  residual  were  tested  against  error,  except  the  combined  regression 
was  compared  against  the  residual  within  environments.  In  the  case  where  both 
heterogeneity  of  regression  and  residual  were  significant,  the  regression  term  was 
tested  against  its  own  residual. 

Stability  parameters 

Stability  parameters  were  calculated  following  Eberhart  and  Russell  (1966). 
Stability  parameter  estimates  such  as  mean,  regression  coefficient,  and  stability 
variance  for  each  genotype  across  environments,  were  computed  to  compare  relative 
stability  of  the  genotypes  for  all  yield  types,  leaf  area,  tiller  dry  weight,  tiller 
density  and  height.  The  stability  variance  was  a  slight  modification  from  Eberhai 
and  Russell's  mean  square  of  deviation  from  regression  (Sd,:).  Stability  variance 
was  measured  by  the  ratio  of  the  deviation  mean  square  (Dev.  MS)  to  the  pooled 
error  mean  square  (MSE/r),  whereas  Sd,2  was  the  difference  between  the  Dev.MS 
and  MSE/r.  A  stable  genotype  is  therefore  defined  as  a  genotype  which  has  an 
above  average  mean  performance,  a  coefficient  of  regression  of  1.0  and  stability 


variance  of  1.0. 
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IV.  RESULTS  AND  DISCUSSION 


SECTION  I 


Environmental  Means 

Genotypic  means  for  each  environment  and  for  all  environments  are  shown  in 
Appendix  Tables  1  to  13  for  all  thirteen  characters.  Means  and  coefficients  of  variation 
(CV)  for  all  genotypes  at  each  of  the  four  locations  and  two  years  are  presented  in 
Tables  7,  8  and  9.  There  are  large  differences  between  locations  and  between  years  for 
yield,  leaf  and  stem  characters.  In  general,  Edmonton  and  Lethbridge  were  recognized 
as  higher  yielding  environments  than  either  Beaverlodge  or  Kinsella  based  on  the 
average  annual  yield  of  the  two  years  (Table  7).  Beaverlodge  is  characterized  by  its 
shorter  and  cooler  summer,  and  Kinsella  is  situated  in  a  drier  part  of  the  province.  The 
field  plot  at  Lethbridge,  though  situated  in  the  semi-arid  southern  part  of  the  province, 
w'as  located  on  irrigated  land.  The  constant  seepage  of  water  from  the  irrigation  ditch 
and  the  infrequent  irrigation  at  the  driest  part  of  the  season  (e.g.  June  1977),  together 
with  its  long  hot  summer,  made  Lethbridge  one  of  the  highest  yielding  areas. 

In  1976  the  highest  mean  spring  and  fall  yields  were  obtained  at  Lethbridge 
and  Edmonton  where  yields  were  approximately  twice  that  harvested  at  Beaverlodge.  In 
1977,  the  best  yield  was  obtained  at  Edmonton  where  yield  was  more  than  two  fold 
that  of  Kinsella.  In  1976  the  fall  yields  were  generally  higher  than  the  spring  yield  at  all 
locations  except  Lethbridge.  A  reverse  pattern  was  displayed  for  1977  due  mainly  to 
the  nitrogen  being  depleted  by  the  vigorous  spring  growth.  The  soil  analysis  conducted 
in  the  fall  (Table  2)  confirmed  that  nitrogen  levels  were  low.  The  low  fall  yield  of  1977 
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Table  7.  Means  and  coefficients  of  variation  (CV)  for  yield  characters  of  smooth  bromegrass  measured  at 
four  locations  and  for  two  years. 
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Table  8.  Means  and  coefficients  of  variation  (CV)  for  leaf  characters  of  smooth  bromegrass  measured  at 
four  locations  and  for  two  years. 
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Table  9.  Means  and  coefficients  of  variation  (CV)  for  tiller  characters  of  smooth  bromegrass  measured  at 
four  locations  and  for  two  years. 
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reflected  that  one  fall  fertilizer  application  per  year  was  inadequate  for  fully  grown 
bromegrass  cut  more  than  once  in  a  season.  Split  fertilizer  applications  would  be 
recommended. 

The  CVs  for  spring,  fall  and  annual  yield  at  Beaverlodge  were  consistently 
higher  than  the  CV's  at  other  locations  in  1976  (Table  7)  indicating  that  the 
measurements  were  subject  to  greater  variation  than  at  other  locations.  This  may  be 
due  in  part  to  the  great  soil  heterogeneity  within  that  area  (Table  1).  The  variation 
became  less  obvious  in  the  1977  data. 

Yields  per  unit  area  were  similar  at  Lethbridge  and  Edmonton,  while 
Beaverlodge  had  the  highest  and  Kinsella  the  lowest  in  1976.  The  same  pattern  holds 
for  1977,  but  the  values  were  smaller.  The  plants  grown  at  Beaverlodge  had  heavier 
leaves  and  stems  (Table  8),  and  more  and  heavier  tillers  (Table  9)  than  at  all  other 
locations  in  both  years.  A  reverse  situation  was  obtained  for  Kinsella  except  that  the 
tiller  density  in  1977  was  high.  All  or  some  of  these  factors  may  have  contributed  to 
the  yield  per  unit  area.  The  interrelationships  among  yield  types  and  between  yield  and 
its  morphological  components  will  be  dealt  with  in  greater  detail  later.  Leaf  number 
and  SLW  were  considered  relatively  stable.  The  leaf  area  per  tiller  was  similar  for  both 
years  at  both  Beaverlodge  and  Edmonton,  but  for  Lethbridge  and  Kinsella,  the  leaf 
area  in  1977  was  only  half  that  of  1976.  Plants  at  Beaverlodge  and  Edmonton  were 
consistently  taller  in  both  years,  and  were  approximately  30%  taller  than  those  of 
Kinsella  and  Lethbridge  in  1977. 

Transformation 

Large  differences  between  the  error  mean  squares  are  shown  for  most  of  the 
characters  measured  in  each  environment  (Tables  10,  11  and  12).  The  relationships 
between  the  means  and  error  mean  squares  are  presented  in  Table  6.  Also  presented  in 
Table  6  are  the  correlation  coefficients  between  the  means  and  error  mean  squares 
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computed  from  the  transformed  data,  showing  that  these  two  statistics  now  are 
independent.  Furthermore,  the  transformation  reduced  the  heterogeneity  of  error  mean 
squares  considerably  for  most  of  the  characters  (Tables  10,  11  and  12).  For  yield  per 
unit  area,  the  difference  between  the  largest  error  mean  square  and  the  smallest  was 
much  greater  on  the  direct  scale  than  on  the  transformed  scale  (Table  10).  For 
example,  the  largest  error  mean  square  for  yield  per  area  on  the  direct  scale  was  32.67 
for  year  1976  at  Beaverlodge,  which  was  18  times  greater  than  the  smallest  error  mean 
square,  1.78  for  1977  in  Edmonton.  However,  on  the  transformed  scale,  the  largest 
error  mean  square  (0.1996)  was  only  5  times  larger  than  the  smallest  mean  square 
(0.04).  For  SLW,  the  largest  error  mean  square  was  82  times  the  smallest  on  the  direct 
scale,  but  2.5  times  the  smallest  on  the  transformed  scale;  and  for  leaf  stem  ratio,  the 
largest  error  mean  square  was  28  times  the  smallest  on  the  direct  scale,  but  5  times  the 
smallest  on  the  transformed  scale.  This  explained  why  transformation  was  applied  to 
SLW  and  leaf  stem  ratio  even  though  there  were  no  apparent  correlations  between 
their  means  and  error  mean  squares  (Table  6).  In  fact,  the  correlation  coefficients 
between  the  means  and  error  mean  squares  for  the  two  characters  were  increased  after 
transformation.  However,  the  smaller  increments  were  compensated  for  by  the 
reduction  in  the  heterogeneity  of  error  mean  squares.  Neither  leaf  number  nor  height 
were  transformed  since  they  showed  either  no  or  only  a  small  correlation  between  the 
mean  and  error  mean  squares  (Table  6).  Also,  their  error  variances  were  sufficiently 
homogeneous  (Tables  11  and  12). 

Analyses  of  variance  were  conducted  on  all  characters  measured  on  both  direct 
and  transformed  scales.  Tests  of  significance  from  computations  using  data  on  the 
direct  scale  were  only  approximations  because  of  the  known  association  between 
evironmental  means  and  error  variances  on  this  scale.  However,  this  is  the  scale 
commonly  used  to  measure  adaptation  parameters  in  agronomic  research,  and  analyses 
using  it  might  be  helpful  for  interpretation  purposes. 
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Table  10.  Error  mean  squares  from  analysis  of  variance  for  individual  environments  for  yield  characters 
measured  on  both  direct  and  transformed  scales. 
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Table  11.  Error  mean  squares  from  analysis  of  variance  for  individual  environments  for  leaf  characters 
measured  on  both  direct  and  transformed  scales. 
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Table  12.  Error  mean  squares  from  analysis  of  variance  for  individual  environments  for  tiller  characters 
measured  on  both  direct  and  transformed  scales. 
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Combined  Analysis  of  Variance  and  Diallel  Analysis 

The  combined  analyses  of  variance  for  yield,  leaf  and  stem  characters  are 
shown  in  Tables  13,  14  and  15,  respectively.  All  the  main  effects  were  highly  significant 
when  tested  by  their  corresponding  errors.  In  general,  the  analyses  show  that  for  both 
scales  the  genotypes  interacted  significantly  with  years  and  with  locations  which  must 
include  climatic,  edaphic  and  management  factors.  The  second  order  interactions, 
genotype  x  location  x  year,  were  significant  for  most  characters,  with  the  exception  of 
stem  dry  weight,  tiller  density  and  leaf  stem  ratio  on  the  direct  scale,  and  leaf  area  and 
leaf  stem  ratio  on  the  transformed  scale.  For  stem  dry  weight  and  tiller  dry  weight,  the 
non-significant  location  by  year  interaction  on  the  transformed  scale  indicated  that  the 
locations  ranked  similarly  from  one  year  to  another.  No  general  pattern  was  obtained 
for  other  characters.  Both  scales  agreed  well  except  in  a  few  instances  which  were 
mainly  the  higher  order  interactions,  where  the  two  scales  were  either  in  disagreement 
or  were  significant  at  different  levels  of  probability. 

As  the  analyses  showed  that  locations,  years  and  their  interactions  were  highly 
significant,  the  data  from  the  four  locations,  two  years,  and  their  interactions  were 
therefore  treated  as  eight  different  environments  in  subsequent  analyses.  Model  I, 
method  2,  of  the  analysis  of  Griffing  (1956)  was  used  to  partition  the  variation  among 
genotypes  into  that  due  to  general  combining  ability  (GCA)  and  specific  combining 
ability  (SC A).  The  interaction  between  genotype  and  environment  was  further  divided 
into  GCA  by  environment  and  SCA  by  environment  interactions  for  all  characters 
(Tables  16,  17  and  18). 

All  GCA  and  SCA  mean  squares  and  their  interactions  with  environments  were 
highly  significant  for  all  characters  on  both  scales.  The  two  scales  closely  agreed  for  all 
cases.  Variation  accounted  for  by  combining  ability  effects,  in  general,  was  appreciably 
higher  than  the  interaction  effects  of  GCA  and  SCA  with  the  environments.  Mean 


Table  13.  Mean  squares  for  the  yield  characters  measured  on  direct  and  transformed  scales  for  the 
seven  parents  and  21  progenies  tested  in  four  locations  and  two  years. 
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Table  14.  Mean  squares  for  the  leaf  characters  measured  on  direct  and  transformed  scales  for  the 
seven  parents  and  21  progenies  tested  in  four  locations  and  two  years. 
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**,  f  significant  at  the  0.05,  0.01,  and  0.001  probability  levels,  respectively. 
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Table  16.  Combining  ability  analysis  for  the  yield  characters  measured  on  direct  and  transformed  scales 
for  the  seven  parents  and  21  progenies  tested  in  eight  environments. 
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Significant  at  the  0.001  probability  level. 


Table  17.  Combining  ability  analysis  for  the  leaf  characters  measured  on  direct  and  transformed  scales 
for  the  seven  parents  and  21  progenies  tested  in  eight  environments. 
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Table  18.  Combining  ability  analysis  for  the  tiller  characters  measured  on  direct  and  transformed  scales 
for  the  seven  parents  and  21  progenies  tested  in  eight  environments. 
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Significant  at  the  0.001  probability  level. 


40 


squares  of  SCA  were  highly  significant  for  all  the  characters.  However,  the  ratio  of 
GCA  to  SCA  showed  that  the  genetic  variability  for  each  character  was  associated  with 
GCA  rather  than  SCA,  indicating  the  importance  of  additive  genetic  effect.  Exceptions 
to  this  were  spring,  fall  and  annual  yield,  where  both  GCA  and  SCA  were  equally 
important  since  the  ratios  were  approximately  2.  The  GCA  to  SCA  ratios  for  the  other 
traits,  on  the  other  hand,  ranged  from  7.8  for  leaf  stem  ratio  to  38.9  for  SLW,  both  on 
the  direct  scale.  These  results  agreed  well  with  those  reported  previously  (Tan  et  al., 
1976b,  1977),  where  results  were  based  on  one  year  and  one  location  only. 

In  contrast  to  the  three  types  of  yield,  the  GCA  to  SCA  ratio  was  high  for 
yield  per  unit  area.  This  may  imply  that  more  rapid  response  would  be  expected 
through  selection  on  yield  per  unit  area  than  selection  on  yield  per  plot. 

The  GCA  effects  of  the  parents  and  the  SCA  effects  of  the  progenies  will  be 
considered  in  detail  in  conjunction  with  stability  parameters  in  Section  II. 

Correlation 

Phenotypic  and  genotypic  correlation  coefficients  among  yield  types  and 
between  yield  and  other  morphological  characters  are  presented  in  Tables  19  and  20, 
respectively.  These  correlation  coefficients  estimated  the  joint  responses  among 
characters  in  this  population.  Phenotypic  correlations  showed  that  the  yield  types  were 
all  significantly  associated.  As  expected,  large  positive  genotypic  correlation  coefficients 
were  obtained  among  all  yield  types. 

On  breaking  down  the  genotypic  correlations  between  characters  into 
additive(GCA)  and  non-additive  (SCA)  genetic  correlation,  the  additive  genetic 
correlation  will  include  all  the  additive  covariance  and  interactions  involving  them,  and 
that  of  non-additive  genetic  correlation  will  include  dominance  and  all  epistatic 
covariances.  The  magnitudes  of  both  additive  and  non-additive  genetic  correlations 
among  yield  types  were  similar  to  or  exceed  the  phenotypic  and  genotypic  correlations 
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Table  19.  Correlation  coefficients  of  phenotypic  (rp) >  genotypic  (r^) 
general  combining  ability  (r  ) ,  and  specific  combining 
ability  (r  )  among  yield  characters  in  smooth  bromegrass. 


Character 

Fall  Yield 

Annual  Yield 

Yield/ area 

Spring  yield 

rP 

0.  78~ 

0.95T 

** 

0.48 

rG 

0.  80 

0.96 

0.  54 

rGCA 

1.09 

1.00 

0.62 

rSCA 

0.  78 

0.96 

0.  72 

Fall  yield 

rP 

0.93T 

0.  67T 

rG 

0.93 

0.  73 

rGCA 

1.02 

1.60 

rSCA 

0.93 

0.76 

Annual  yield 

rP 

0.  58T 

rG 

0.66 

rGCA 

0.92 

rSCA 

0.  79 

§  Statistical 

test  is 

provided  for  r^ 

only. 

**,  f  significant  at  the  0.01,  and  0.001  probability  levels, 
respectively,  where  n=28. 
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(Table  19).  The  additive  genetic  correlations  were  considerably  greater  than  the 
non-additive  genetic  correlations  in  most  cases,  except  between  spring  yield  and  yield 
per  unit  area. 

Very  few  of  the  morphological  characters  measured  were  associated  with  yield 
(Table  20).  Height  was  positively  while  leaf  stem  ratio  was  negatively  correlated  with 
spring,  fall  and  annual  yield.  Tiller  density  was  significantly  correlated  phenotypically 
with  all  yield  types  except  fall  yield.  Significant,  negative  correlation  coefficients  were 
noted  for  SLW  and  spring  and  annual  yield.  More  morphological  characters  showed 
significant  phenotypic  correlations  with  yield  per  area  than  with  other  yield  types. 

These  morphological  characters  included  leaf  area,  leaf,  stem  and  tiller  dry  weight,  tiller 
density  and  height.  Genotypic  correlation  coefficients  were  of  similar  or  slightly  higher 
magnitudes  than  the  phenotypic  correlation  coefficients.  Further  partition  of  the 
genotypic  response  showed  that  the  close  genotypic  associations  between  the  height  and 
the  spring  and  annual  yields  were  due  mainly  to  the  non-additive  genetic  effects  and 
little,  if  any,  to  the  additive  genetic  effects.  The  close  associations  between  tiller  density 
and  the  types  of  yield  were  mainly  due  to  additive  genetic  variance,  while  the 
association  due  to  non-additive  genetic  variation  were  negative,  small  and  negligible. 

The  large  additive  genetic  correlation  coefficients  indicated  that  in  this  population 
selection  towards  high  yield  would  be  in  conjunction  with  increased  tiller  density.  For 
high  yield  per  area  selection  through  large  leaf  area,  greater  number  of  tiller  and  per 
tiller  dry  weight  should  also  be  effective. 

The  number  of  tillers  per  unit  ground  area  (or  tiller  density,  in  no/dm")  and 
tiller  dry  weight  are  considered  as  the  two  major  components  of  forage  yield.  The 
relative  importance  of  these  two  components  in  their  contribution  to  forage  yield  have 
been  reported  previously  by  various  workers.  In  the  present  study  and  the  previous 
reports  (Tan  et  al.,  1977a,  b)  tiller  density  showed  a  highly  positive  correlation  with 
forage  yield,  but  the  association  between  tiller  dry  weight  and  forage  yield  was  not 
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significant.  The  results  are  in  close  agreement  with  those  in  Lolium  species  (Lazenby 
and  Rogers,  1962;  Thomson  et  al.,  1973),  and  suggest  that  tiller  density  is  of  greater 
importance  in  determining  yield  than  tiller  dry  weight.  Conversely,  different  results  on 
the  associations  between  tiller  density,  tiller  size  and  yield  have  been  reported  in  a 
number  of  forage  species,  including  orchardgrass  (Knight,  1970),  perennial  ryegrass 
(Troughton,  1965),  and  bromegrass  (Walton,  1976). 


SECTION  II 


Joint  Regression  Analysis 

Tables  21,  22  and  23  show  the  further  partition  of  the  GE  interactions  using 
the  joint  regression  analysis  on  both  scales.  In  this  analysis  regressions  of  genotype 
means  from  environments  (i.e.  locations  and  years)  were  compared  with  environmental 
means  derived  from  the  whole  population.  Such  analyses  made  possible  the 
determination  for  each  character  how  much  of  the  GE  interaction  could  be  accounted 
for  by  the  heterogeneity  of  the  regressions,  that  is  by  the  differences  between  slopes  of 
the  regression  lines,  and  how  much  was  residual  and  therefore  unpredictable.  The  two 
components  were  tested  against  the  error  terms  derived  from  analysis  of  variance. 
Regardless  of  which  character  was  being  examined  a  significant  part  of  the  GE 
interaction  could  always  be  demonstrated  as  being  due  to  the  differences  between  the 
slopes  of  the  regression  lines.  In  all  instances  the  unpredictable  residual  deviations  were 
also  significant  on  both  scales. 

When  both  the  mean  square  of  heterogeneity  of  regression  and  the  mean 
square  of  residual  were  significant  when  tested  against  the  error  term,  the  validity  of 
any  predictions  would  depend  on  the  relative  magnitudes  of  the  two  mean  squares 
(Perkins  and  Jinks,  1968).  If  the  mean  square  of  heterogeneity  of  regression  was 
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Table  22.  Mean  squares  from  the  joint  regression  analyses  for  the  leaf  characters  measured  on  direct  and 
transformed  scales  for  28  genotypes  tested  in  eight  environments,  with  environmental  means 
derived  from  all  genotypes. 
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at  the  0. 05  probability  level  for  the  heterogeneity  of  regression  when  compared 


Table  23.  Mean  squares  from  the  joint  regression  analyses  for  the  tiller  characters  measured  on  direct 

and  transformed  scales  for  28  genotypes  tested  in  eight  environments,  with  environmental  means 
derived  from  all  genotypes. 
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significant  when  compared  with  that  of  residual,  the  prediction  of  GE  interactions 
based  on  the  linear  regression  would  still  have  considerable  practical  value.  Both  yield 
per  area  and  fall  yield  showed  significant  regression  when  compared  with  the  residual 
on  both  scales  (Table  21).  Variations  due  to  heterogeneity  among  regressions  for  leaf 
area  and  leaf  dry  weight  were  both  significant  on  the  direct  scale  but  not  on  the 
transformed  scale  (Table  22).  Mean  squares  for  heterogeneity  of  regressions  were 
significant  for  stem  dry  weight,  tiller  dry  weight,  tiller  density  and  leaf  stem  ratio  on 
the  direct  scale,  but  only  stem  dry  weight  and  tiller  density  were  significant  on  the 
transformed  scale. 

Thus,  the  linear  model  would  retain  considerable  predictive  value  for  the 
genotypes  concerned  for  some  of  the  characters  such  as  tiller  density,  stem  dry  weight, 
yield  per  area  and  fall  yield  as  most  of  the  interactions  were  accounted  for  by  the 
heterogeneity  of  regressions.  According  to  Perkins  and  Jinks  (1968),  even  if  the  mean 
squares  due  to  heterogeneity  among  regressions  was  not  significant  when  tested  against 
residual  mean  square;  this  did  not  rule  out  the  possibility  that  the  regressions  of  each 
genotype  mean  on  to  the  environmental  mean  for  some  of  the  genotypes  taken 
individually  may  be  highly  significant  when  tested  against  their  residual  mean  squares. 
Therefore,  for  those  particular  genotypes  reliable  predictions  can  still  be  made. 

Attention  to  the  problems  that  can  arise  from  assessing  an  environment  by  the 
mean  of  the  genotypes  grown  in  it  has  been  drawn  by  Freeman  and  Perkins  (1971). 
Their  main  criticism  of  the  technique  was  that  the  genotype  means  contribute  to,  and 
hence  were  not  statistically  independent  of,  the  environmental  means  on  which  they 
were  regressed.  Yates  and  Cochran  (1938)  originally  justified  this  procedure  on  the 
ground  that,  because  the  joint  regression  sum  of  squares  equalled  the  environmental 
sum  of  squares,  the  GE  interaction  sum  of  squares  could  then  be  partitioned 
orthogonally  into  two  items,  one  measuring  differences  between  the  slopes  of  the  fitted 
lines,  the  other  measuring  the  residual  deviations  about  these  lines.  Freeman(  1973) 
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agreed  that  this  procedure  was  perfectly  valid  provided  inferences  were  drawn  only 
about  the  sample  of  genotypes  and  environments  used  in  the  experiment,  that  is  the 
model  was  fixed. 

Assuming  a  fixed  model,  the  joint  regression  analysis  presented  in  Tables  21, 

22  and  23  are  considered  statistically  valid.  The  present  trial  includes  seven  parents  and 
21  progenies  in  all  environments,  which  coincides  with  one  of  the  experimental  designs 
outlined  by  Freeman  and  Perkins  (1971),  who  approached  these  problems  by  providing 
independent  estimates  of  the  environment  by  the  mean  of  the  parental  genotypes.  In 
this  study  the  method  of  Freeman  and  Perkins(  1971)  was  followed  and  the  results  of 
the  analysis  are  presented  in  Tables  24,  25  and  26  for  yield,  leaf,  and  tiller  characters, 
respectively,  on  both  scales. 

Most  of  the  terms  in  Table  24  are  significant  at  the  0.1  per  cent  level  when 
compared  with  the  error.  The  mean  squares  for  heterogeneity  of  regressions  were 
significantly  greater  than  their  residuals  for  yield  per  area,  spring  yield  and  fall  yield  on 
the  direct  scale,  but  were  significant  only  for  spring  and  fall  yield  on  the  transformed 
scale.  For  each  character,  the  combined  regression  slope  was  not  significantly  different 
from  unity,  but  the  residual  about  this  regression  was  not  negligible.  Mean  squares  for 
heterogeneity  among  regressions  were  significantly  greater  than  the  error  for  all  leaf 
characters  except  SLW  on  both  scales  (Table  25).  When  compared  with  their  residuals, 
mean  squares  of  the  heterogeneity  of  regression  were  significant  for  leaf  area  and  leaf 
dry  weight  on  the  direct  scale,  but  not  on  the  transformed  scale. 

In  Table  26,  height  was  the  only  character  which  did  not  show  a  significant 
linear  response  on  the  direct  scale  when  regressions  were  compared  with  either  error  or 
residual  terms.  The  mean  squares  of  heterogeneity  among  regressions  for  tiller  density 
remained  significant  while  leaf  stem  ratio  was  not  significant  after  the  scale  had  been 
transformed. 

The  effect  of  transformation  has  become  obvious  in  the  joint  regression 


Table  24.  Mean  squares  from  the  joint  regression  analyses  for  the  yield  characters  measured  on  direct  and 
transformed  scales  for  21  genotypes  tested  in  eight  environments,  with  environmental  means 
derived  from  the  parental  genotypes. 
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Table  25.  Mean  squares  from  the  joint  regression  analyses  for  the  leaf  characters  measured  on  the 
and  transformed  scales  for  21  genotypes  tested  in  eight  environments,  with  environmental 
derived  from  the  parental  genotypes. 
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analysis.  Test  of  significance  among  means  discussed  in  Section  I  were  similar  on  both 
scales.  For  heterogeneity  among  regressions,  the  tests  produced  different  results.  This 
suggested  that  tests  of  significance  on  the  direct  scale  were  accurate  among  means,  but 
were  not  accurate  among  regressions.  Cases  where  drastic  changes  of  variance  ratios 
when  different  scales  are  used  were  not  unusual.  Transformation  may  be  used  to  reduce 
both  interaction  effects  and  heterogeneity  of  variances  (Bartlett,  1947).  However,  if  the 
GE  interaction  could  be  removed  by  using  any  of  the  transformation  encountered  in 
this  study  (or  any  other  monotonic  transformation),  the  interaction  will  likely  be  of  the 
kind  that  cannot  be  exploited  by  selection  (Eagles  et  al.,  1977). 

The  use  of  different  measures  of  environment  also  created  some  differences  in 
the  significance  of  the  linear  responses.  For  example,  most  of  the  interaction  for  spring 
yield  was  accounted  for  by  heterogeneity  of  regressions  on  both  scales  when  the 
environmental  means  were  measured  by  the  mean  of  the  parental  genotypes,  but  was 
not  significant  on  either  scale  when  the  overall  mean  of  the  28  genotypes  was  taken 
into  the  measure  of  environment.  These  conflicting  results  from  the  present  study 
cautions  the  choice  of  measure  of  environment.  It  would  affect  the  interpretation  of 
results  from  regression  analysis.  The  differences  which  arose  from  the  two  measures  of 
environments  would  likely  be  the  result  of  the  use  of  too  few  parental  genotypes  as  a 
measure  of  environment,  as  postulated  by  Perkins  and  Jinks  (1973). 

Stability  Parameters 

Stability  parameter  estimates  for  each  genotype  across  environments  were 
computed  to  compare  relative  stability  of  the  genotypes.  A  stable  genotype  was  defined 
as  a  genotype  which  had  a  regression  coefficient  of  1 .0  and  stability  variance  around 
1.0.  Of  all  traits,  yield  characters  such  as  spring,  fall  and  annual  yields  are  of  major 
agronomic  interest  and  are  the  most  important  criteria  on  which  stability  of  field  crops 
is  being  considered.  Yield  per  area  and  other  morphological  characters,  such  as  leaf 
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area,  tiller  weight  and  density,  and  height  were  also  included,  to  show  that  stability  of 
genotypes  may  vary  with  traits,  that  is,  the  genotypes  which  were  rated  as  stable  for 
one  character  may  be  unstable  for  other  characters. 

The  primary  stability  parameter  for  plant  cultivars,  according  to  Finlay  and 
Wilkinson  (1963),  is  the  regression  coefficient  of  the  means  of  a  cultivar  upon  the 
environmental  indexes.  This  statistic  measures  whether  a  given  cultivar  produces 
relatively  better  or  worse  than  all  cultivars  in  the  environments  with  low-  and 
high-productivity  potentials.  Herein,  b— 1.0  will  be  designated  as  average  stability  (see 
Finlay  and  Wilkinson,  1963),  values  above  unity  will  be  classed  as  less  stable  than 
average  (below  average  stability)  and  values  below,  as  more  stable  than  average  (above 
average  stability).  Since  the  heterogeneity  among  regressions  was  not  a  very  large 
proportion  of  the  GE  interaction  for  most  characters,  the  partition  of  the  pooled 
deviation  into  stability  variances  for  every  genotype  appeared  very  important.  It  is 
therefore  necessary  to  consider  the  two  stability  parameters  together  with  mean 
performances  of  each  genotype,  mainly  following  Eberhart  and  Russell's  (1966) 
definition  of  stability,  viz.,  6,  =  1.0  and  Sd,2=0. 

The  stability  parameters  of  each  genotype,  including  parents  and  progenies,  for 
spring,  fall  and  annual  yields  are  given  in  Tables  27,  28  and  29  on  both  scales.  Graphic 
summaries  are  also  given  in  Figures  1,  2  and  3  for  spring,  fall  and  annual  yields, 
respectively,  on  direct  scale  only  since  those  of  the  transformed  scale  closely  followed 
similar  patterns.  The  graphs  show  the  relationships  between  yield  (spring,  fall  and 
annual)  and  linear  response  of  genotypes,  and  are  considered  useful  in  selecting  stable 
genotypes.  The  vertical  line  is  the  grand  mean,  whereas  the  horizontal  line  is  the 
average  slope  ( 6=  1 .0).  All  stability  variances  significantly  greater  than  the  error 
variance  (P<0.05)  are  shown  by  symbol  7^'.  For  clarity,  a  code  number  instead  of 
pedigree,  as  given  in  Table  4,  is  marked  for  each  genotype. 

The  spring  yield  varied  from  1655  to  2335  grams  per  plot.  The  regression 
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Table  27.  Means,  regression  coefficients,  and  stability  variances  for 
spring  yield  measured  on  both  direct  and  transformed  scales 
on  28  bromegrass  genotypes. 


Genotypes 

Direct  Scale 

Transformed 

Scale 

Mean 

Regression 

coefficient 

S  tability 
variance 

Regression 

coefficient 

Stability 

variance 

Parents 

(g) 

1 

2198 

1.19^ 

1.30 

1.18 

1.18 

2 

2107 

0.82 

5. 02 1 

0.  81 

4 . 54t 

3 

1970 

0.98 

2.  77* 

1.02 

2.97** 

4 

1655 

0.95 

6. 91f 

1.03 

8. 06f 

5 

1726 

1.09 

3.  77f 

1.15 

3. 50** 

6 

2316 

1.18 

3. 71** 

1.13 

2.99** 

7 

2071 

1.  21 

10. 99 1 

1.19 

11. 20t 

Progenies 

12 

2317 

0.  83 

2.00 

0.80 

1.80 

13 

2269 

0.92 

1.07 

0.92 

1.07 

14 

1984 

1.15 

0.95 

1. 17f 

0.  89 

15 

1969 

1.03 

4 . 62f 

1.03 

4. 52f 

16 

2288 

1.07 

1.81 

1.05 

1.69 

17 

2205 

0.  88 

1.96 

0.  88 

1.66 

23 

1882 

0.  72f 

1.  78 

0.  75  + 

1.  78 

24 

2102 

0.94 

0.  70 

0.95 

0.66 

25 

2141 

1.01 

1.02 

1.00 

1.02 

26 

2269 

0.  86 

0.  63 

0.  84  f 

0.54 

27 

2246 

0.91 

1.  79 

0.  87 

1.68 

34 

2335 

1.12 

2.08 

1.10 

1.68 

35 

1948 

0.94 

0.69 

0.94 

0.  78 

36 

2119 

0.97 

0.  65 

0.95 

0.68 

37 

2135 

0.77 

4. 60f 

0.  77 

3. 90t 

45 

2199 

1.25 

2.68* 

1.  25 

2.39* 

46 

1864 

1.01 

1.07 

1.06 

1.06 

47 

19  82 

1.01 

3.19** 

1.05 

2.87** 

56 

1954 

1.07 

5. 20f 

1.06 

4.38+ 

57 

2268 

1.05 

0.56 

1.03 

0.  58 

67 

2326 

1.06 

2.90** 

1.02 

2.32* 

Grand  mean 

2102. 

1 

LSD  (.05) 

237 

denotes  regression  coefficients  that  were  higher  or  lower  than 
1.0  at  the  0.05  probability  level. 


f  denote  deviations  from  regression  greater  than  the  pooled  error 
in  analysis  of  variance  at  the  0.05,  0.01,  and  0*001  probability 
levels,  respectively. 
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Fig.  1.  The  relationship  between  spring  yield  (g/plot)  and  stability 
of  28  bromegrass  genotypes. 

Stability  variance  significant  (P  <  0.05); 

+  Stability  variance  not  significant. 
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Table  28.  Means,  regression  coefficients,  and  stability  variances  for 
fall  yield  measured  on  both  direct  and  transformed  scales 
on  28  bromegrass  genotypes. 


Genotypes 

Direct 

Scale 

Transformed  Scale 

Mean 

Regression 

coefficient 

S  tability 
variance 

Regression 
coef  f  icient 

Stability 

variance 

Parents 

(g) 

1 

1546 

1.06 

2.04 

1.08 

3. 56** 

2 

1643 

1.11 

3.  86  + 

1.09 

2.39* 

3 

1468 

0.  76 

3.04** 

0.78 

2.98** 

4 

1336 

1.13 

5.73+ 

1.13 

6.28+ 

5 

1386 

0.  62 

8.05  + 

0.64 

8.12  + 

6 

1803 

1.10 

1.07 

1.02 

0.89 

7 

1475 

1.18 

5.65  + 

1.25 

3.92  + 

Progenies 

12 

1752 

0.99 

1.58 

0.92 

1.47 

13 

1718 

0.99 

0.58 

0.92 

0.40 

14 

1431 

0.  98 

1.  32 

1.03 

1.08 

15 

1605 

0.  78 

5.  51  + 

0.  77 

4.60  + 

16 

1691 

0.96 

0.  81 

0.94 

0.  88 

17 

1614 

1. 32  + 

1.52 

1.27  + 

1.42 

23 

1381 

0.  75+ 

0.  73 

0. 79f 

0.58 

24 

1581 

1.04 

0.31 

1.02 

0.35 

25 

1655 

0.  90 

1.65 

0.91 

1.28 

26 

1666 

0.  98 

1.14 

0.97 

0.79 

27 

1576 

1.19 

1.71 

1.17 

1.47 

34 

1727 

1.17 

2.87** 

1.13 

2.17* 

35 

1594 

0.91 

0.47 

0.  88* 

0.32 

36 

1648 

0.96 

1.56 

0.95 

1.39 

37 

1656 

1.10 

2.04 

1.03 

2.09 

45 

1440 

0.  98 

0.  33 

1.04 

0.  64 

46 

1159 

0.  87  + 

0.  61 

1.03 

1.  26 

47 

12  89 

1.  07 

1.  50 

1.15 

1.70 

56 

1484 

0.  82 

2.17* 

0.82 

2.09 

57 

1663 

1.14 

2.35* 

1.11 

1.96 

67 

1610 

1. 16 

2.60* 

1.15 

2.08 

Grand  mean 

1557. 

5 

LSD  (.05) 

205 

f 

deno  tes 

regression 

coefficients 

that  were  higher 

or  lower  than 

1.0  at  the  0.05  probability  level. 

,  t  denote  deviations  from  regression  greater  than  the  pooled  error 
in  analysis  of  variance  at  the  0.05,  0.01,  and  0.001  probability 
levels,  respectively. 
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Fig-  2.  The  relationship  between  fall  yield  (g/plot)  and  stability 
of  28  bromegrass  genotypes. 

)K  Stability  variance  significant  (P  <  0.05); 

+  Stability  variance  not  significant. 
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Table  29.  Means,  regression  coefficients,  and  stability  variances  for 
annual  yield  measured  on  both  direct  and  transformed  scales 
on  28  bromegrass  genotypes. 


Genotypes 

Direct  Scale 

Tr ansformed 

[  Scale 

Mean 

Regression 

coefficient 

Stability 

variance 

Regression 

coefficient 

Stability 
var iance 

Parents 

(g) 

1 

3744 

1.16 

2.09 

1.15 

2.25* 

2 

3750 

0.  84 

4.44  + 

0.85 

3.  87  + 

3 

3439 

0.  82 

3.93+ 

0.90 

4.51  + 

4 

2992 

1.14 

6.75+ 

1.21 

9.36+ 

5 

3112 

0.92 

8.99  + 

0.96 

8.37+ 

6 

4119 

1.12 

2.73* 

1.07 

2.09 

7 

3547 

1.  27 

9.56+ 

1.25 

8.  82  + 

Progenies 

12 

4070 

0.  81 

1.92 

0.  77 

2.04 

13 

3988 

0.93 

0.  78 

0.  88 

0.  72 

14 

3415 

1.12 

1.01 

1.16 

0.99 

15 

3574 

0.  81 

6 . 54  + 

0.82 

5.  74  + 

16 

3980 

1.01 

1.71 

0.99 

1.  55 

17 

3819 

1.09 

3.73** 

1.06 

2. 97** 

23 

3263 

0.  70+ 

1.40 

0.  76  + 

1.  35 

24 

3684 

0.  99 

0.  68 

0.99 

0.  60 

25 

3796 

0.93 

1.  59 

0.92 

1.49 

26 

3936 

0.  85+ 

0.  36 

0.  82 

0.40 

27 

3822 

1.  02 

2.54* 

0.99 

2.61* 

34 

4063 

1.26 

2.19* 

1.18 

1.68 

35 

3542 

0.  88 

0.63 

0.90 

0.  62 

36 

3767 

1.07 

0.64 

1.05 

0.  59 

37 

3791 

0.90 

4.  87  + 

0.88 

4.20+ 

45 

3639 

1.15 

1.93 

1.16 

1.91 

46 

3023 

0.94 

0.  55 

1.01 

0.  74 

47 

3271 

1. 12 

2.40* 

1.18 

1.94 

56 

3438 

0.92 

4 . 53+ 

0.96 

3.45** 

57 

3932 

1.04 

1.33 

1.02 

1.26 

67 

3937 

1.16 

3. 13** 

1.09 

2.51** 

Grand  mean 

3659. 5 

LSD  (.05) 

400 

+  denotes  regression  coefficients  that  were  higher  or  lower  than 

1.0  at  the  0.05  probability  level. 


*,  **,  f  denote  deviations  from  regression  greater  than  the  pooled  error 
in  analysis  of  variance  at  the  0.05,  0.01,  and  0.001  probability 
levels,  respectively. 
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Fig.  3.  The  relationship  between  annual  yield  (g/plot)  and  stability 
of  28  bromegrass  genotypes. 

*  Stability  variance  significant  (P  <  0.05); 

+  Stability  variance  not  significant. 
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coefficients  ranged  from  0.72  to  1.25  (Table  27  and  Fig.  1).  Four  genotypes,  namely  1, 

14,  23  and  26  had  regression  coefficients  significantly  different  from  unity  on  one  scale 
or  both  (Table  27).  Six  out  of  the  seven  parental  genotypes  and  six  progeny  genotypes 
had  significant  deviation  mean  squares.  Joint  consideration  of  parameters  indicated  that 
genotypes  1,  12,  13,  16,  17,  24,  25,  27,  34,  36,  37,  45,  57,  6  and  67  combined  high  yield 
and  average  linear  response.  Genotypes  37,  45,  6  and  67,  however,  recorded  significant 
deviations.  Genotype  26  had  high  yield  and  above  average  linear  response(6<1.0) 

(Table  27)  and  genotype  1,  on  the  other  hand,  had  below  average  linear  response 

(6>1.0). 

Fall  yield  varied  from  1159  (genotype  46)  to  1803  (genotype  6)  grams  per  plot 
and  the  regression  coefficient  ranged  from  0.65  (genotype  5)  to  1.32  (genotype  17) 
(Tables  28  and  Fig.  2).  Seventeen  genotypes  recorded  higher  yield  than  the  mean  of  all 
genotypes.  These  genotypes  were:  6,  12,  34,  13,  16,  57,  37,  26,  25,  2,  36,  17,  67,  15,  35, 
24  and  27.  All  these  genotypes  except  17  and  35  had  average  linear  response  (6=1). 
Deviation  mean  squares  were  significant  for  five  genotypes,  namely  34,  57,  2,  67  and 

15.  This  suggested  general  adaptability  of  most  of  the  high  yielding  genotypes  for  fall 
yield. 

The  annual  yield  varied  from  2992  (genotype  4)  to  4119  (genotype  6)  grams  per 
plot,  and  the  regression  coefficients  ranged  from  0.7  (genotype  23)  to  1.27  (genotype  7) 
(Table  29  and  Fig.  3).  Sixteen  genotypes  had  higher  yield  than  the  mean  ol  all 
genotypes.  Among  these  genotypes  only  one  (i.e.  26)  had  regression  coefficient  which 
was  significantly  smaller  than  unity,  and  seven,  namely  6,  34,  67,  17,  27,  37  and  2,  had 
significant  deviation  mean  squares. 

Joint  consideration  of  the  three  yield  types  (spring,  fall  and  annual)  and  the 
stability  parameters  indicated  that  five  genotypes  had  consistently  combined  high  yields, 
average  linear  responses  and  non-significant  deviations.  These  genotypes  were;  12,  13, 

16,  25,  and  36.  Genotype  23,  on  the  other  hand,  had  shown  maximum  stability 
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combined  with  low  yields.  Genotype  26  had  combined  high  yield  and  above  average 
response  for  spring  and  annual  yield,  and  average  linear  response  for  fall  yield.  This 
genotype  may  be  preferred  to  genotype  23  where  a  highly  stable  genotype  is  desired. 

Correlation  coefficients  between  regression  coefficients  and  mean  yields  of  each 
genotype  was  0.08  for  spring  yield.  Positive  but  non-significant  correlations  were 
obtained  between  regression  coefficients  and  mean  fall  yields  (r=0.33)  and  between 
regression  coeficients  and  annual  yields  (r^O.H).  Therefore,  it  seems  possible  to  select 
genotypes  which  have  above  average  stability  and  high  yields,  such  as  genotype  26. 
Thess  results  are  in  contrast  to  those  reported  by  Eberhart  and  Russell  (1966)  in  maize, 
and  by  Fatunla  and  Frey  (1976)  and  Eagles  et  al.  (1977)  in  oats,  where  they  obtained 
significant  correlations  between  mean  yields  and  regressions. 

Tables  30  to  34  are  presented  to  show  that  a  phenotype  which  is  classified  as 
stable  for  one  character  may  be  unstable  for  other  characters.  For  yield  per  area  (Table 
30)  none  of  the  genotypes  which  had  high  yield  per  area  were  close  to  average  stability. 
Table  31  shows  that  among  those  genotypes  (1,  16,  25,  35,  36,  57,  and  67)  which  had 
large  leaf  areas  and  were  rated  as  having  average  stability,  three  genotypes  (16,  25  and 
36)  also  had  average  stability  in  yields  (spring,  fall  and  annual).  For  tiller  dry  weight 
(Table  32),  genotypes  1,  12,  13,  14,  34,  35,  and  36  had  average  stability,  but  among 
them  only  genotypes  13,  35  and  36  combined  higher  tiller  weight  than  their  overall 
mean.  Genotype  37  with  average  stability  for  fall  yield  was  the  only  genotype  which 
had  average  stability  for  high  tiller  density  (Table  33).  Genotypes  13,  15,  16,  34,  35,  36, 
56  and  67  combined  both  average  stability  and  above  average  plant  height  (Table  34). 

The  data  indicate  that  it  is  rather  difficult  to  select  a  genotype  which  had 
average  stability  for  yields  and  will  also  be  of  average  stability  for  other  morphological 
characters,  such  as  high  tiller  density,  heavy  tiller  weight,  large  leaf  area,  and  plant 
height.  Genotype  36  probably  combined  average  stability  and  high  means  for  the  most 
characters,  which  included  yields  (spring,  fall  and  annual),  leaf  area,  tiller  dry  weight 


63 


Table  30.  Means,  regression  coefficients,  and  stability  variances  for 

yield  per  area  measured  on  both  direct  and  transformed  scales 
on  28  bromegrass  genotypes. 


Genotypes 

Direct  Scale 

Transformed 

Scale 

Mean 

Regression 

coefficient 

Stability 

variance 

Regression 

coefficient 

Stability 

variance 

Parents 

(§/ dm^) 

1 

37.  7 

0.99 

2.79* 

0.97 

5.06  + 

2 

39.4 

1.  23f 

6.  28+ 

1.08 

5.06  + 

3 

33.1 

0.  69f 

17.41+ 

0.  83 

10.12+ 

4 

23.4 

0.  52+ 

6.65  + 

0.76 

16.63+ 

5 

35.1 

0.  82f 

6.74  + 

0.80 

8.67  + 

6 

40.4 

1.09 

2.75* 

0.95 

2.16* 

7 

25.1 

0.  53f 

14.59+ 

0.88 

23. 86+ 

Progenies 

12 

43.  3 

1.48+ 

12. 22+ 

1.13 

6.50+ 

13 

41.  7 

1.19  + 

6.73+ 

1.00 

13.73+ 

14 

34.7 

1.07 

1.06 

1.10  + 

1.44 

15 

37.0 

1.  04 

2.06 

0.97 

2.16* 

16 

41.  5 

1.33+ 

7.50  + 

1.07 

3.61** 

17 

30.  5 

0.  89+ 

0.40 

1.02 

0.72 

23 

37.  8 

1.11 

5.07+ 

0.97 

5.  78+ 

24 

35.  5 

1.07 

2.90** 

1.11 

5.78+ 

25 

36.  8 

1.19  + 

6.02  + 

1.07 

4.33  + 

26 

40.0 

1.23 

11.03+ 

1.08 

5.06  + 

27 

32.0 

1.  02 

9.67  + 

1.07 

9.39  + 

34 

34. 8 

1.02 

5.04  + 

1.12 

9.39  + 

35 

36.3 

1.15+ 

1.  54 

1.08 

1.44 

36 

38.9 

1.  35  + 

9.08+ 

1.13 

4.33  + 

37 

30.6 

0.  84+ 

2.  72* 

0.99 

6.50  + 

45 

33.  6 

0.97 

3.12** 

1.02 

4.33+ 

46 

25.  8 

0.  71  + 

5.18+ 

1.01 

10.84+ 

47 

23.  7 

0.58+ 

2.07 

0.  89 

7.23  + 

56 

37.5 

1.11 

6.  60  + 

1.00 

+- 

oo 

• 

m 

57 

31.0 

0.  79  + 

7.25  + 

0.95 

6.  50+ 

67 

33.9 

0.94 

10. 53+ 

0.94 

9.39  + 

Grand  mean 

34.7 

LSD  (.05) 

5. 1 

f  denotes  regression  coefficients  that  were  higher  or  lower  than 

1.0  at  the  0.05  probability  level. 


*,  **,  t  denote  deviations  from  regression  greater  than  the  pooled  error 
in  analysis  of  variance  at  the  0.05,  0.01,  and  0.001  probability 
levels,  respectively. 
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Table  31.  Means,  regression  coefficients,  and  stability  variances  for 
leaf  area  measured  on  both  direct  and  transformed  scales  on 
28  bromegrass  genotypes. 


Genotypes 

Direct  Scale 

Transformed  Scale 

Regression 
Mean  coefficient 

Stability 

variance 

Regression 

coefficient 

S  tability 
variance 

Parents 

,  2  v 

(cm  ) 

1 

89.9 

0.91 

2.06 

0.92 

2.03 

2 

76.  8 

0.  84 

1.29 

0.91 

1.42 

3 

69.  4 

0.  75 

6.  80+ 

0.  84 

8.17+ 

4 

70.  6 

0.  79+ 

0.44 

0.91 

0.50 

5 

114.4 

1.13 

3.75+ 

0.99 

2.69* 

6 

118.  7 

1.  58* 

1. 89 

1.39* 

1.37 

7 

57.  2 

0.  71 

2.40* 

0.90 

3.43** 

Progenies 

12 

86.1 

0.  88 

1.28 

0.90 

1.31 

13 

94.9 

0.  8lf 

0.47 

0.  76* 

0.47 

14 

73.4 

0.  80f 

0.  52 

0.  89 

0.  58 

15 

118.2 

1.32 

4.43+ 

1.16 

3.03** 

16 

107.4 

1.  20 

1.  77 

1.08 

1.45 

17 

74.3 

0.93 

0.  25 

1.04 

0.18 

23 

74.  6 

0.68* 

0.  58 

0.  74* 

0.  71 

24 

69.1 

0.  83 

1.30 

0.93 

1.31 

25 

99.9 

1.15 

1.  07 

1.10 

1.05 

26 

107.0 

1.22f 

0.  85 

1.10 

0.  71 

27 

69.  7 

0.  80 

0.  82 

0.91 

0.  89 

34 

83.  8 

0.90 

1.  53 

0.96 

1.95 

35 

101.  7 

1.08 

1.  53 

1.01 

1.26 

36 

100.  2 

1.16 

1.99 

1.11 

1.63 

37 

67.5 

0.  87 

0.  39 

1.02 

0.34 

45 

112. 2 

1.  42* 

1.18 

1.  27* 

0.  84 

46 

87.  7 

0.94 

2.58* 

0.96 

2.08 

47 

69.3 

0.  88 

1.  72 

1.02 

2.11* 

56 

111.  3 

1.  29 

2.58* 

1.17 

1. 89 

57 

92.7 

1.07 

1.96 

1.03 

1.55 

67 

90.0 

1.06 

0.  80 

1.  05 

0.  73 

Grand  mean 

88.9 

LSD  (.05) 

7.  9 

*  denotes  regression  coefficients  that  were  higher  or  lower  than 

1.0  at  the  0.05  probability  level. 

*,  **,  f  denote  deviations  from  regression  greater  than  the  pooled  error 
in  analysis  of  variance  at  the  0.05,  0.01,  and  0.001  probability 
levels,  respectively. 
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Table  32.  Means,  regression  coefficients,  and  stability  variances  for 
tiller  dry  weight  measured  on  both  direct  and  transformed 
scales  on  28  bromegrass  genotypes. 


Direct  Scale  Transformed  Scale 


Genotypes 

Mean 

Regression 

coefficient 

Stability 

variance 

Regression 

coefficient 

Stability 

variance 

(g) 

Parents 

1 

2.66 

0.91 

1.45 

0.  88 

2.00 

2 

2.12 

0.89 

2.04 

1.11 

3.45** 

3 

2.15 

0.  55  + 

5.70+ 

0.  73 

5.50  + 

4 

2.09 

0.  78+ 

0.52 

0.98 

0.  95 

5 

4.05 

1.64 

6.95  + 

1.20 

4.  27  + 

6 

3.  69 

1.60  + 

3.46** 

1.22 

1.77 

7 

1.  87 

0.  53  + 

2.83** 

0.  83 

4.70  + 

Progenies 

12 

2.45 

1.03 

0.  89 

1.07 

1.17 

13 

2.  70 

0.99 

1.07 

0.96 

1.31 

14 

2.12 

0.95 

0.  69 

1.18 

1.00 

15 

4.13 

1.  51 1 

1.69 

1.11 

1.00 

16 

3.  28 

1.26  + 

2.00 

1.04 

1.40 

17 

2.23 

0.  82  + 

0.  79 

0.96 

1.07 

23 

1.  99 

0.77  + 

1.45 

0.93 

3.12** 

24 

2.00 

0.64  + 

0.  27 

0.  83+ 

0.50 

25 

3.35 

1.14  + 

0.34 

0.96 

0.32 

26 

3.00 

1.  22+ 

0.48 

1.11 

0.  54 

27 

1.97 

0.  85  + 

0.  38 

1.06 

0.50 

34 

2.42 

0.93 

1.00 

1.05 

1.15 

35 

3.53 

1.17 

1.59 

0.96 

1.32 

36 

3.02 

1.18 

1.34 

1.07 

1.12 

37 

1.  94 

0.  80+ 

0.20 

1.  05 

0.  37 

45 

3.  80 

1.21+ 

1.17 

0.96 

0.  70 

46 

2.43 

0.66  + 

2.76* 

0.  80 

2.90** 

47 

2.  00 

0.  74  + 

1.90 

0.98 

. 2.87** 

56 

3.  58 

1.32  + 

0.44 

1.07 

0.37 

57 

3.14 

0.88 

2.28* 

0.  81 

1  .87 

67 

2.84 

1.06 

2. 62* 

1.05 

2.40* 

Grand  mean  2.73 
LSD  (.05)  0.32 


+  denotes  regression  coefficients  that  were  higher  or  lower  than 

1.0  at  the  0.05  probability  level. 

*,  **,  t  denote  deviations  from  regression  greater  than  the  pooled  error 
in  analysis  of  variance  at  the  0.05,  0.01,  and  0.001  probability 
levels,  respectively. 
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Table  33.  Means,  regression  coefficients,  and  stability  variances  for 

tiller  density  measured  on  both  direct  and  transformed  scales 
on  28  bromegrass  genotypes. 


Genotypes 

Direct  Scale 

Tr ansf ormed 

1  Scale 

Mean 

Regression 

coefficient 

Stability 

variance 

Regression 

coefficient 

S  tability 
variance 

(no /dm2) 

Parents 

1 

28.1 

1 . 46f 

8.08 

1.43t 

7.97  + 

2 

28.3 

1.51 

11. 38+ 

1.42 

10. 69+ 

3 

27.1 

1.11 

7.65  + 

1.07 

6.99  + 

4 

17.  8 

0.44t 

11.74+ 

0.51 

14.21+ 

5 

16.1 

0.56 

23.46+ 

0.  68 

30.40+ 

6 

21.9 

0.85 

12.74+ 

0.  88 

11.38+ 

7 

21.  7 

0.  87 

9.02  + 

0.  89 

9.71  + 

Progenies 

12 

28.1 

1.  77t 

5.  25  + 

1. 60t 

3. 58** 

13 

25.2 

1.20 

4.99  + 

1.15 

5.08+ 

14 

25.2 

1.15 

5.49  + 

1.12 

5.31  + 

15 

16.2 

0.  7 It 

2.40* 

0.84 

3.12** 

16 

23.2 

0.93 

4.11  + 

0.91 

3.  69** 

17 

21.  6 

1.06 

2.30* 

1.07 

2.71* 

23 

31.9 

1.45 

10.98+ 

1.  29 

7.97  + 

24 

26.  2 

1.20 

14.47+ 

1.15 

14.04+ 

25 

19.1 

0.  79 1 

1.40 

0.  86 

1.73 

26 

25.  2 

1.09 

2.56* 

1.05 

2.19* 

27 

25.1 

1.11 

12.63+ 

1.03 

10.57+ 

34 

22.6 

1.07 

3.14** 

1.09 

3.35** 

35 

18.0 

0.  81 

4.93  + 

0.91 

5.  66  + 

36 

25.4 

1.20 

7.13  + 

1.15 

6.82  + 

37 

25.  8 

0.96 

1.90 

0.90 

1.67 

45 

16.3 

0.  83 

2.63* 

0.95 

3. 69** 

46 

18.  6 

0.92 

2.49* 

0.99 

2.60* 

47 

19.4 

0.75 

11.34+ 

0.  79 

12.48+ 

56 

19.3 

0.  80 

3.12** 

0.87 

4.04  + 

57 

17.6 

0.69t 

1.50 

0.  77f 

2.08 

67 

23. 2 

0. 67f 

1.84 

0.68t 

1.  79 

Grand  mean 

22.5 

LSD  (.05) 

2.5 

t  denotes  regression  coefficients  that  were  higher  or  lower  than 

1.0  at  the  0.05  probability  level. 


*,  ** ,  T  denote  deviations  from  regression  greater  than  the  pooled  error 
in  analysis  of  variance  at  the  0.05,  0.01,  and  0.001  probability 
levels,  respectively. 
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Table  34.  Means,  regression  coefficients,  and  stability  variances 
for  height  of  28  broraegrass  genotypes. 


Genotypes 

Direct  Scale 

Mean 

Regression 

coefficient 

Stability 

variance 

(cm) 

Parents 

1 

111 

0.  97 

2.42* 

2 

101 

1.  25 

2.86** 

3 

105 

1.04 

8  .  79 1 

4 

97 

0.  76 

4. 33t 

5 

114 

0.92 

4 . 18f 

6 

127 

0.99 

4 . 08f 

7 

97 

1.06 

5. 44t 

Progenies 

12 

107 

0.97 

2.56* 

13 

113 

0.94 

1.48 

14 

101 

1.04 

2.59* 

15 

116 

1.09 

1.04 

16 

118 

1.00 

1.88 

17 

105 

1.00 

3.09** 

23 

99 

0.96 

2.43* 

24 

101 

0.90 

0.50 

25 

113 

1.181= 

0.  70 

26 

113 

1.07 

3.72** 

27 

100 

1.10 

2.24* 

34 

108 

0.92 

0.99 

35 

114 

1.01 

1.33 

36 

114 

1.04 

1.81 

37 

102 

1.01 

6.  57  + 

45 

113 

0.98 

2.50* 

46 

96 

0.  88 

5.48+ 

47 

95 

0.  88 

2.73* 

56 

114 

0.95 

1.10 

57 

107 

1.01 

2.45* 

67 

113 

1.04 

1.65 

Grand  mean 

108.1 

LSD  (.05) 

4 

^  denotes  regression  coefficients  that  were  higher  or  lower 

than  1.0  at  the  0.05  probability  level. 

*,  **,  t  denote  deviations  from  regression  greater  than  the  pooled 
error  in  analysis  of  variance  at  the  0.05,  0.01,  and  0.001 
probability  levels,  respectively. 
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and  plant  height.  It  seems  therefore  reasonable  to  conclude  that  different  genotvpes 
could  be  classified  as  stable  for  different  sets  of  characters  even  though  all  the 
characters  were  measured  from  the  same  set  ot  genotypes. 

There  were  also  differences  between  the  two  scales  used  in  the  tests  of 
significance  for  both  regression  coefficients  and  stability  variances.  In  no  cases  will  a 
genotype  be  considered  stable  unless  the  parameters  of  both  scales  agree  to  certain 
extent. 

Relationships  Between  Stability  Parameters  and  Combining  Abilities  for  Yields 

The  diallel  design  of  the  present  experiment  permits  the  partition  of  regression 
coefficients  and  deviation  mean  squares  into  GCA  and  SCA.  Therefore,  the  combining 
ability  of  a  genotype  for  both  parameters  and  yield  can  be  jointly  examined. 

Combining  ability  analyses  for  the  means,  linear  regression  and  deviation  mean  squares 
for  the  yield  types  are  given  in  Table  35.  Only  the  three  yield  types,  namely  spring,  fall 
and  annual  yields,  are  presented  in  the  joint  consideration  between  combining  ability  of 
mean  yield  and  combining  ability  for  the  stability  parameters,  i.e.  regression  coetlicient 
and  deviation  mean  square,  with  respect  to  practical  selection  in  breeding. 

Variation  in  the  expression  of  mean  yield  due  to  GCA,  SCA  and  their 
interactions  with  environments  was  significant  when  compared  with  the  pooled  error 
(Table  35).  The  relative  magnitudes  of  GCA  to  SCA  indicated  both  GCA  and  SCA 
were  equally  important  for  mean  yields.  For  regression  coefficient,  the  variation  due  to 
GCA  was  significantly  greater  than  that  of  SCA.  This  applied  to  all  three  > ield  t>pes. 
The  GCA  variance  of  deviation  mean  square  was,  on  the  other  hand,  not  significantlv 
greater  than  that  of  SCA.  Hence,  the  differences  in  linear  response  of  genotypes  to 
environments  were  primarily  due  to  GCA,  while  both  GCA  and  SCA  were  equally 
important  in  the  expression  of  deviations.  Data  on  both  scales,  direct  and  transformed 
(viz.  power  and  logarithm),  closely  agreed  for  the  levels  of  significance.  Furthermore, 


Table  35.  Mean  squares  for  combining  ability  of  stability  parameters:  means  (X),  regression  coefficients 
(b)  and  deviation  mean  squares. 
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the  direct  scale  was  the  scale  commonly  used  to  measure  adaptation  parameters  in 
agronomic  research.  Therefore,  only  the  direct  scale  will  be  presented  in  subsequent 
discussion. 

The  average  GCA  effects  for  spring  yield  ranged  from  -115  to  71  (Fig.  4).  Four 
parents,  1,  2,  6  and  7,  had  GCA  effects  greater  than  zero;  parent  6  being  the  best 
combiner  for  spring  yield.  The  estimates  of  GCA  for  the  regression  coefficient  varied 
from  -0.121  to  0.059  (Fig.  4).  Parent  5  had  the  highest  GCA  effects  for  regression 
coefficient  of  spring  yield;  while  parent  7  had  near  zero  GCA  effect.  The  lowest  GCA 
effect  was  obtained  for  parent  2,  followed  by  parent  3.  Figure  4  shows  the  association 
between  GCA  effects  of  mean  yield  and  GCA  effects  of  regression  coefficient.  Parents  6 
and  1,  the  best  general  combiners  for  mean  spring  yield  performance,  had  positive 
GCA  effects  for  regression  coefficient.  Parent  7  seems  promising  since  it  had 
moderately  high  GCA  for  spring  yield  and  near  zero  GCA  for  regression.  Thus, 
according  to  the  average  GCA  effect  presented  in  Figure  4,  parent  7  would  transmit  to 
its  progeny  high  spring  yield  and  average  stability.  Parent  2  transmitted  above  average 
stability  to  its  progeny  and  had  positive  GCA  estimates  for  spring  yield. 

Figure  5  examines  the  association  of  GCA  effects  between  spring  yield  and 
regression  with  respect  to  individual  environments.  The  figure  shows  apparent  GCA  x 
environment  interactions  by  the  different  ranking  of  the  environments  within  each 
genotype.  It  revealed  that  the  environmental  mean  of  a  genotype  as  shown  in  Figure  4 
was  not  as  informative  as  that  shown  by  individual  environments  when  genotype  x 
environment  interaction  was  prevalent.  For  example,  parent  7  was  among  the  highest 
combiners  for  spring  yield  when  average  GCA  for  all  environments  was  considered 
(Fig.  4),  but  when  individual  environments  were  examined  (Fig.  5),  it  had  positive 
GCA  effects  only  at  four  of  the  eight  environments.  Parent  2  transmitted  to  its  progeny 
high  spring  yield  of  above  average  stability  in  five  of  the  eight  environments.  Parent  1 
had  positive  GCA  effects  for  spring  yield  at  the  highest  number  of  environments  -  all 
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Fig.  4.  The  relationship  between  mean  GCA  effects  for  spring  yield 
and  regression  coefficient  for  all  environments. 
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Fig.  5.  The  relationship  between  GCA  effects  for  spring  yield  and 
regression  coefficient  for  eight  individual  environments. 
(B : Beaverlodge ;  K:Kinsella;  L:Lethbridge;  ErEdmonton) 
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eight,  and  a  near  average  stability.  It  was  followed  by  parent  6,  which  had  positive 
GCA  effects  for  six  environments,  and  positive  GCA  effects  for  regression  coefficients. 

With  the  exception  of  parent  1  which  combined  the  highest  spring  yield  and 
average  stability,  no  generalization  can  be  made  for  the  rest  of  the  parents.  It  appears 
rather  difficult  to  obtain  a  parent  which  combined  well  at  all  environments. 

For  fall  yield,  parents  1,  2  and  6  transmitted  average  stability  to  their  progenies 
and  had  the  highest  GCA  effects  for  mean  fall  yield  (Fig.  6).  Parents  5  and  7  had, 
respectively,  the  lowest  and  the  highest  GCA  for  regression  coefficients.  Figure  7 
provides  information  on  the  GCA  effects  of  fall  yield  and  regression  coefficient  for 
each  environment.  Parent  1  and  6  combined  equally  well  for  fall  yield  and  average 
stability  at  all  environments.  For  parent  2,  there  was  a  greater  variation  among 
environments  than  for  the  aforementioned  parents.  With  the  exception  of  one 
environment,  i.e.  Lethbridge  1976,  parent  2  transmitted  above  average  fall  yield  and 

average  stability  to  its  progeny. 

Considering  annual  hay  production  (Fig.  8),  parents  1  and  6  again  had  higher 
GCA  effects  than  the  rest  of  the  parents;  the  former  had  near  average  stability  while 
the  latter  had  rather  high  GCA  effect  for  regression.  Both  parents  showed  high 
consistencies  over  all  environments  (Fig.  9).  Parent  2  transmitted  above  average 
stability  as  well  as  high  annual  yield  to  its  progeny  at  five  of  the  eight  environments 

studied. 

joint  consideration  of  the  GCA  effects  of  the  three  yield  characters  indicated 
that  parent  I,  a  random  selection  from  ‘Magna'  which  is  one  of  the  leading  cultivars  in 
Western  Canada,  seems  to  have  the  highest  potential  among  the  seven  clones.  As 
shown  in  Figures  1,  2  and  3,  the  parent  I  itself  was  neither  the  best  yielder  nor  stable. 
Yet  as  indicated  by  the  GCA  effects,  progenies  of  this  clone  will  likely  have  high  yield 
and  average  stability,  and  will  show  consistent  performance  over  environments.  The 
next  good  yield  combiner,  parent  6  had  positive  GCA  effects  for  linear  response.  This 
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Fie.  6.  The  relationship  between  mean  GCA  effects  for  fall  yield 
and  regression  coefficient  for  all  environments. 
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Fig.  7.  The  relationship  between  GCA  effects  for  fall  yield  and 

regression  coefficient  for  eight  individual  environments. 
(B : Beaverlodge ;  K:Kinsella;  L :Lethbridge ;  E:Edmonton) 
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Fig.  8.  The  relationship  between  mean  GCA  effects  for  annual  yield 
and  regression  coefficient  for  all  environments. 
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The  relationship  between  GCA  effects  for  annual  yield  and 
regression  coefficient  for  eight  individual  environmen  ts . 
(B : Beaverlodge ;  K:Kinsella;  L :Lethbridge ;  E:Edmonton) 
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parent  may  have  some  potential  value  in  areas  which  are  fertile  and  properly  managed. 
Parent  2  demonstrated  the  highest  GCA  effects  for  stability  of  all  parents;  however,  it 
did  not  combine  equally  high  yield  in  such  favourable  environments,  as  Edmonton  in 
1977,  Lethbridge  in  both  1976  and  1977,  as  in  other  less  favourable  environments 
(Figure  9).  Parent  3  transmitted  above  average  stability  but  rather  low  yield  to  its 
progenies  in  many  of  the  environments  tested. 

Figures  10,  11  and  12  show  the  association  between  GCA  effects  of  mean 
yields  (spring,  fall  and  annual)  and  the  GCA  effects  for  mean  squares  of  deviations 
from  regression.  These  figures  again  confirm  that  parents  1  and  6  combined  high  yield 
and  small  deviation  mean  squares  with  high  consistency  over  the  range  of 
environments.  Parents  2  and  3  were  the  next  best  in  their  abilities  to  transmit  small 
deviation  from  regression  to  their  progenies,  but  were  inconsistent  in  yielding 
performances  in  varying  environments.  Parents  5  and  7  had  the  highest  mean  squares 
of  deviation,  hence  were  highly  unstable  for  transmission  of  yielding  performances  to 
their  progenies. 

The  GCA  effects  for  mean  yields  and  the  GCA  effects  for  linear  response  were 
not  significantly  correlated  in  any  of  the  yield  types  (r=-0.26  for  spring;  r  0.11  for  fall, 
and  r-- 0.35  for  annual  yield).  These  suggest  a  lack  ot  association  between  these 
attributes  and  the  possibility  of  breeding  cultivars  which  combined  high  yield 
performance  and  stability. 

Since  both  GCA  and  SCA  were  equally  important  in  determining  yields,  the 
associations  between  SCA  effects  of  mean  yield  and  SCA  effects  of  regression 
coefficients  were  examined  in  details  in  Figures  13  to  18  inclusive,  for  the  three  yield 
types  on  environmental  means  (Figs.  13,  15  and  17)  and  on  individual  environments 

(Figs.  14,  16  and  18). 

Figure  13  shows  that  genotype  34  combined  the  highest  mean  SCA  effect  ot 
spring  yield  and  regression  coefficient,  and  was  followed  by  45,  67,  16  and  13. 
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Fig.  10.  The  relationship  between  GCA  effects  for  spring  yield  and 
deviation  mean  square  for  eight  individual  environments. 
(B : Beaverlodge;  K:Kinsella;  L:Lethbridge;  E:Edmonton) 
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Fie.  11.  The  relationship  between  GCA  effects  for  fall  yield  and 
deviation  mean  square  for  eight  individual  environments. 
(B : Beaver lodge;  L: Lethbridge;  K:Kinsella;  E:Edmonton) 
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Fig.  12.  The  relationship  between  GCA  effects  for  annual  yield  and 
deviation  mean  square  for  eight  individual  environments. 
(B : Beaverlodge;  KrKinsella;  LrLethbridge;  E:Edmonton) 
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Fig.  13.  The  relationship  between  mean  SCA  effects  for  spring  yield 
and  regression  coefficient  for  all  environments. 

*  P  <  0.05. 
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Fig  14.  The  relationship  between  SCA  effects  for  spring  yield  and 
regression  coefficient  for  eight  individual  environments. 
(B : Beaverlodge ;  K:Kinsella;  L rLethbr idge ;  ErEdmonton) 
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Fig.  15.  The  relationship  between  mean  SCA  effects  for  fall  yield 
and  regression  coefficient  for  all  environments. 

***  p  <  0.001. 
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Fig.  16.  The  relationship  between  SCA  effects  for  fall  yield  and 

regression  coefficient  for  eight  individual  environments. 
(B : Beaver  lodge ;  K:Kinsella;  L : Lethbridge ;  ErEdmonton) 
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Fie  18  The  relationship  between  SCA  effects  for  annual  yield  and 
regression  coefficient  for  eight  individual  environments. 
(B : Be aver lodge ;  K:Kinsella;  L:Lethbridge;  E:Edmonton) 
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Genotypes  12  and  26,  on  the  other  hand,  had  low  mean  SCA  estimates  for  regression 
and  above  average  SCA  for  spring  yield.  Genotype  57  combined  high  yield  and  average 
stability.  Figure  14  shows  that  all  of  these  genotypes  showed  certain  consistencies  over 
environments.  Genotype  34  again  showed  the  highest  mean  SCA  effect  for  both  fall 
yield  and  regression  (Fig.  15)  and  had  the  highest  consistency  in  all  eight  environments 
(Fig.  16).  Genotypes  16,  67  and  12  had  above  average  SCA  for  fall  yield  and  average 
stability.  None  of  the  genotypes  combined  consistently  high  fall  yield  and  above 
average  stability.  Genotype  34  again  had  remarkably  high  mean  SCA  estimates  for 
both  annual  yield  and  regression  (Fig.  17)  and  displayed  the  highest  consistency  in  all 
environments  (Fig.  18).  Genotypes  57,  16,  24  and  13  had  higher  than  average  SCA  for 
annual  yield  and  near  zero  mean  SCA  for  regression.  Genotype  12  combined  high  yield 
and  above  average  stability  in  most  of  the  environments  (Fig.  18).  Genotypes  23,  46 
and  47  had  combined  above  average  stability  and  consistently  low  yield.  There  were 
highly  significant  correlation  coefficients  between  mean  SCA  effects  of  regression  and 
spring  (r=0.54,  0.05  level),  fall  (r=0.78,  0.001  level),  and  annual  (r=0.68,  0.001  level) 
yields.  Therefore,  genotypes  which  had  high  SCA  estimates  for  yield  tended  to  be 
relatively  unstable,  and  those  which  had  low  SCA  for  yields  tended  to  be  stable. 

The  results  obtained  from  the  association  of  SCA  effects  between  yields  and 
regression  were  in  general  agreement  with  those  obtained  from  the  graphs  showing  the 
relationship  between  mean  yields  and  regression  (Figs.  1,  2  and  3),  although  the  latter 
pairs  showed  no  significant  correlations.  Genotypes  12,  13,  16  and  24,  which  were 
identified  earlier  as  potentially  good  genotypes  were  again  found  to  have  desirable 
potential  by  combining  ability  analysis.  Genotypes  25  and  36  appeared  more  desirable 
in  the  mean-regression  association  than  in  the  combining  ability  analysis.  Genotype  34 
had  its  potential  in  favourable  environments;  unfortunately,  because  its  stability 
variances  were  significant  for  both  fall  and  annual  yields,  certain  limits  would  be 
imposed  on  the  use  of  this  otherwise  superior  genotype.  Genotype  23  was,  at  the  other 


extreme,  stable  but  had  a  yield  performance  which  was  far  too  low  to  be  useful  in 
practical  breeding. 
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V.  GENERAL  DISCUSSION 


In  the  combined  analysis  of  variance,  the  main  effects  of  years,  locations  and 
genotypes  were  highly  significant  for  yield,  leaf  and  tiller  characters  in  smooth 
bromegrass.  The  analysis  also  showed  that  the  genotypes  interacted  significantly  with 
locations  and  years.  Since  genotypes  and  environmental  interactions  were  highly 
significant,  no  immediate  generalization  could  be  made  on  the  relative  performance  of 
these  genotypes  over  even  a  restricted  range  of  environments.  Since  the  analysis  of 
variance  gave  no  further  useful  explanation  of  the  GE  interaction,  the  data  were  futher 
analysed  by  the  joint  regression  methods  as  proposed  by  Perkins  and  Jinks  (1968). 

The  joint  regression  analysis  by  ascribed  a  significant  part  of  the  GE 
interaction  to  the  heterogeneity  of  the  regression  lines  in  every  character.  The  regression 
technique  turned  complex  GE  interactions  into  a  series  of  linear  and  thus  predictable 
responses.  However,  the  residual  component  was  also  significant  in  most  instances,  so 
there  was  a  certain  amount  of  unpredictable  and  unaccountable  variation. 

The  model,  proposed  by  Finlay  and  Wilkinson  (1963)  and  further  elaborated 
by  Eberhart  and  Russell  (1966),  provided  a  method  of  screening  individual  genotypes. 
They  suggested  that  regression  coefficients  and  deviation  mean  squares  be  considered 
with  mean  performance  as  criteria  for  stability.  However,  it  is  difficult  to  decide  the 
relative  weight  to  be  attached  to  these  three  parameters  while  selecting  materials  in  a 
practical  breeding  program.  Some  workers  (Eberhart  and  Russell,  1966,  1969;  Busch, 
Hammond  and  Frohberg,  1976)  have  suggested  that  the  deviation  mean  square  is  a 
more  important  stability  parameter  than  linear  regression. 

In  the  present  investigation,  both  the  heterogeneity  among  regressions  and  the 
residuals  were  highly  significant  for  most  of  the  characters.  The  relative  magnitudes  of 
mean  squares  suggested  that  in  most  cases  where  regressions  were  not  significantly 
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greater  than  their  residuals,  both  linear  regression  and  deviation  mean  square  should  be 
considered  for  each  genotype. 

Results  from  the  two  methods  of  assessment  of  environments,  namely,  by  the 
means  of  all  genotypes  and  of  parental  genotypes,  led  some  different  conclusions  for 
some  of  the  characters  measured  in  the  present  studies.  Fripp  and  Caten  (1971)  found 
little  differences  between  the  use  of  a  control  genotype  and  the  use  of  the  mean  of  a 
large  number  of  the  test  genotypes  to  estimate  environments.  Perkins  and  Jinks  (  1973) 
also  obtained  similar  values  of  significance  between  environmental  means  measured  by 
a  large  number  of  genotypes  and  those  derived  from  other  closely  related  sets  of 
genotypes.  Regressions  on  mean  derived  from  only  a  few  independent  genotypes  might 
be  the  cause  of  the  problems  of  interpretation  (Perkins  and  Jinks,  1973). 

Diallel  analysis  indicated  that  both  GCA  and  SCA  were  involved  and  of  equal 
importance  in  the  inheritance  of  yields,  including  spring,  fall  and  annual  yields. 
Interactions  (GCA  x  E  and  SCA  x  E)  also  accounted  for  significant  variation,  thereby 
indicating  the  importance  of  multi-environmental  testing.  For  regression  estimates, 
mean  squares  due  to  GCA  were  about  3.5  (spring  and  annual  yields)  to  5.6  (fall  yield) 
times  as  great  as  those  due  to  SCA;  while  mean  squares  due  to  GCA  and  SCA  for 
deviations  did  not  differ  significantly.  A  predominance  of  additive  effects  in  the 
inheritance  of  linear  responses  was  recorded  by  Eberhart  and  Russell  (1966,  1969)  and 
Patanothai  and  Atkins  (1974).  Both  additive  and  non-additive  effects  were  involved  in 
inheritance  of  deviations  from  regression  (Eberhart  and  Russell,  1969;  Dhillon  and 
Singh,  1977).  Busch  et  al.  (1976),  however,  indicated  that  deviation  from  regression  was 
simply  inherited  and  could  be  predicted  from  the  parents.  In  contrast,  Patanothai  and 
Atkins  (1974)  concluded  that  the  inheritance  of  deviations  from  regression  was 
complex. 

The  presence  of  a  substantial  proportion  of  variability  due  to  the  additive 
genetic  component  in  the  inheritance  of  linear  response  suggests  that  it  should  be 
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possible  to  exploit  this  fraction  of  variability  in  developing  high  yielding  stable 
cultivars.  However,  the  differences  among  the  genotypes  for  mean  yield  were  ascribable 
to  both  additive  and  non-additive  genetic  variance.  Therefore,  a  more  complicated 
approach,  such  as  recurrent  selections  involving  multi-location  and  multi-year  testing, 
seems  necessary. 

None  of  the  parents  had  general  adaptability  and  high  mean  performance  for 
all  three  yield  types.  Parent  1,  however,  had  high  and  consistent  mean  GCA  effects  for 
mean  yields,  near  zero  GCA  effects  for  linear  response  and  low  GCA  effects  for 
deviations  from  regression,  thus  transmitting  both  high  yields  and  average  stability.  It 
may  be  the  most  suited  parent  for  hybridization  program.  However,  since  SC  A  was 
equally  as  important  as  GCA  in  determining  yield,  not  all  the  progenies  of  parent  1 
were  equally  high  in  yield.  This  seemed  obvious  because  only  half  of  the  crosses 
(genotypes  12,  13  and  16)  involving  parent  1  as  one  of  their  parents  had  above  average 
yield  as  well  as  average  stability.  One  other  good  combiner,  parent  6,  had  positive 
GCA  effects  for  linear  response  which  suggested  that  its  progeny  can  do  well  in  fertile 
and  well-managed  environments,  but  not  as  well  in  unfavourable  environments. 

The  results  presented  here  illustrated  that  regression  analyses  are  powerful  tools 
in  the  analysis  of  GE  interactions.  In  cases  where  deviation  mean  squares  were  as 
important  as  the  heterogeneity  among  regressions,  further  partition  of  the  pooled 
deviation  into  individual  genotypes  proved  useful  in  discriminating  genotypes  for 
stability  parameters.  The  technique  has  been  used  to  analyse  yield  in  a  number  of  other 
experiments  involving  different  grass  population  and  species  from  the  smooth 
bromegrass  reported  here.  In  addition,  stability  of  characters  other  than  yield  was  also 
evaluated,  which  is  rarely  done  in  other  crops  except  peas  (Snoad  and  Arthur,  1974). 
When  yield  is  considered  alone,  two  successive  harvests  and  the  total  annual  yield  were 
included.  Only  those  genotypes  which  showed  some  degree  of  consistency  for  all  three 
yield  types  were  selected.  The  use  of  diallel  cross  design  in  the  present  experiment 


93 


permits  the  joint  investigation  of  combining  ability  and  the  stability  parameters. 

With  the  diverse  climatic  conditions  in  Alberta,  the  data  showed  that  to  obtain 
genotypes  which  will  be  consistently  high  in  yield  and  highly  stable  for  the  whole 
province  would  be  difficult.  One  of  the  alternatives  would  be  to  sub-divide  the  whole 
province  into  zones.  Under  such  division,  more  genotypes  are  expected  to  show  specific 
adaptation  as  well  as  high  mean  performance.  Because  of  the  few  locations  used  in  this 
study,  no  immediate  suggestion  could  be  made  for  the  possible  zoning  for  the  province. 
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APPENDICES 


Appendix  Table  1.  Mean  spring  yield  (g/plot)  of  smooth  bromegrass  grown  at  eight  environments  in  Alberta. 
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Appendix  Table  2.  Mean  fall  yield  (g/plot)  of  smooth  bromegrass  grown  at  eight  environments  in  Alberta. 
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Appendix  Table  3.  Mean  annual  yield  (g/plot)  of  smooth  bromegrass  grown  at  eight  environments  in  Alberta. 
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Appendix  Table  4.  Mean  yield  per  area  (g/dm  )  of  smooth  bromegrass  grown  at  eight  environments  in  Alberta. 
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Appendix  Table  5.  Mean  leaf  area  (cm  /tiller)  of  smooth  bromegrass  grown  at  eight  environments  in  Alberta. 
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Appendix  Table  6.  Mean  leaf  number  (no/ tiller)  of  smooth  bromegrass  grown  at  eight  environments  in 
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Appendix  Table  7.  Mean  leaf  dry  weight  (g/tiller)  of  smooth  bromegrass  grown  at  eight  environments  in  Alberta. 
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Appendix  Table  9.  Mean  stem  dry  weight  (g / tiller)  of  smooth  bromegrass  grown  at  eight  environments 


109 


o 

•H 


a 

>v  c 
■u  ctf 


O 

e 

<u 

o 


a) 

s 


,  r-  ^  i  -  -  ^^NiNCO^Osjmillcsli^HffiOHrHOvO  m  CO 

rHr'.r'rvOcNOinocNr-'CovDcovovorN.j-vojavav'Ci-vooa'ivoa'Nvocn  cncsi 


r^roa\<rcr>iooooo^H 


CN 


I— lr— IHmrOt— ICNCNr- I  CO  CN  r- It— IrlMNHrlNNHrOr- IHCNCNCM 


CM  O 


fi 

O 

•u 

G 

o 

£ 

x> 

w 


r^ 

f^N 

CTV 


vO 

CTV 


.98 

.95 

r^ 

co 

.41 

.46 

.46 

On 

st 

o 

CM 

.06 

00 

<3- 

ON 

CN 

rH 

st 

oo 

in 

.51 

.51 

o 

m 

st 

CN 

rH 

St 

r^ 

m 

74 

v£> 

St 

rs 

st 

m 

vO 

CN 

OV 

35 

75 

12 

IT 

CO 

O 

m 

St 

T - 1 

rH 

i — I 

rH 

CN 

CN 

r — 1 

rH 

CN 

i — 1 

CO 

CN 

rH 

rH 

rH 

CN 

CN 

rH 

rH 

CM 

CN 

i — 1 

CN 

1 - 1 

rH 

CM 

CN 

CN 

CN 

o 

CO 

m 

00 

CN 

rs 

st 

CM 

vO 

iH 

in 

o 

m 

00 

i — 1 

CM 

St 

ao 

»H 

CTn 

rs 

is. 

00 

CO 

st 

co 

LT) 

• 

rH 

LO 

00 

co 

CM 

o 

CN 

st 

CM 

00 

rs 

CN 

VO 

o 

CO 

co 

st 

rs 

rs 

st 

OV 

rH 

CN 

fs 

CO 

in 

st 

CM 

CM 

rH 

co 

CM 

CN 

CM 

CN 

CO 

CM 

CM 

rH 

rH 

CO 

CN 

rH 

CM 

CO 

CM 

rH 

co 

l — 1 

CM 

CO 

CM 

CM 

CN 

o 

a) 

bC 

T3 


r^. 

ov 

rH 


u 

rO 

X 

4J 

0) 

hJ 


VO 

av 


fO  i — I  LO  , — |  CO  vf  CM  I — |  LP|  VO 

n.  i — I  i — I  i— I  lO  CO  i — |  lT)  i — I 

1  I  i“t  rH  rH  CM  CM  t — I  i — i  i — I  i — I 


CNCMvOrHCOO'iOOrHOCN 

r^Mt-iOvooor'.r^vocovo 

' — I  ' — IrHi— (CMCNrHi— li— |i— I 


a\inHo>NvoHiOfOov 

'd-OvcMCTNCNOCOiH'd-O 


CM)  i — I  i— I  O  rH  CN  1 — I  I — |  I — |  CM 


r^cOiHocOMOcoooo 

OfOM>f\OiOHiON<r 

CO  CM|  i — I  i— I  1 — I  CM  CM  ' — I  I — |  CM) 


vococococor^ocrv  co  <f 

CO  O  fC|  1C|  r- I  CM  CM  VO  CO  CO 
•••••••*  •• 

i — I  i — I  CM  H  H  CM  i — I  i — I  i — IO 


vocoomcoovcoovo  oo  is 
i — I  in  i — I  i — I  iO  CO  IM  cO  o  O 

CMHCOCMHCMCMN  CMJO 


I — i 


r^. 

OM 


H I  CO  CM  VO  in  CO 
VO  OV  CM  OM  CO  N 

i — I  O  i — I  O  i — I  i— I 


O  O  CO  IO  M  M 
M  CM  M-  CO  H  N 

O' — I  ' — I  O  CN  t— H 


o  o  vo  cn  m  av 

H  CM  Oc  Oc  M-  CO 
'  I  1 - 1  O  '  I  1“H  O 


H  vO  M  00  CN  CM 
H  OC  IO  00  H  H 


i — I  ' — I  rH  O  CM  i — I 


co  oc  im  ac  oo  in 

CO  VO  00  CN  CO  CN 

O  rH  i — |  i — |  , — |  O 


0) 

CO 

c 

•H 


vO 

r^. 

OV 


OOOOOONf 

lOiHMCOIMM- 

i-H  H  H  rH  CM  CN 


CN  LO  CO  O  NT  vO 
CO  in  CO  H  CO  rH 

' — I  i — I  i — I  H  CN  CN 


<f  O  co  o  CTv 

IO  H  n  o  M  H 

1 — I  rH  i — I  CN  ' — I  ' — I 


m  r-^  vo  co 

CO  sf  O  CN  M 


' — I  CN  i — I  i —  CM  i — I 


CO  o  vO  CN  co  o 

CO  CN  CN  CO  P's  CO 

H  CM  CN  H  rH  O 


CO 


0) 

n 


a) 

bd 

x) 

o 


I"'- 

f''. 

OV 


rH  vO  vO  P"^  r-H 
O  O  o 

CM  CM  CM  CM  M 


N  H  M"  CO  CO 

LO  OM  Mp  VO  rH 

CO  rH  CN  CM  CN 


v£)  IM  ov  CO  CO 

CO  CN  O  CN  O 

Vt  CO  CN  CN  CM 


O  <t  <t  CO 
CO  o  O  VO  vO 


CO  CO  CN  CM  CO 


icuo  CO  o\  O 
ON  O  vO  vf  ^ 


CN  CN  CO  CN  rH 


in  on  vf  co  cn 

io  cm  o  ~ct 

CO  CO  CO  CN  O 


>-! 

<D 

> 

cfl 

a) 

PQ 


vO 

I--. 

OM 


o  m  vo  on  m  in 

vo  o  co  vo 

co  co  CN  cn  m  in 


CO  in  H  M  H  O 
rH  VO  CO  O  m  00 

cn  co  co  co  m 


CO  <f  CO  CO  CO  CO 

o  co  m  in  co  o 


O  CM  CM  M  Vf  co 


ON  H  vf  O  CN  CN 
i — I  vO  CO  CO  W  (m 


CO  -<p  Sj"  CM  vt  CM 


im  vo  co  vt  cm  in 

in  co  co  cm  p^  in 

CN  O’  CO  CO  CO  O 


<v 

Cu 


o 

1=1 

0) 

o 


HCMCOsJiOnOImcN 


CO  Vf  IT)  vO 


MCO-JIOnONMIOvDMiOvOMvOMM 

HCMCMCMCNCMCOCOCOCOM-vtM-iOiONO 


|X 


in 

o 


Q 

cn 

r4 


Appendix  Table  10.  Mean  tiller  dry  weight  (g/tiller)  of  smooth  bromegrass  grown  at  eight  environments 

in  Alberta. 
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Appendix  Table  11.  Mean  tiller  density  (no /dm  )  of  smooth  bromegrass  grown  at  eight  environments  in 

Alberta. 
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Appendix  Table  12.  Mean  leaf  stem  ratio  of  smooth  bromegrass  grown  at  eight  environments  in  Alberta. 
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Appendix  Table  13.  Mean  plant  height  (cm)  of  smooth  bromegrass  grown  at  eight  environments  in  Alberta. 
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